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ABSTEACT 

A fundamental study of viscoseals having a rarefied gas as the 
sealant has been conducted. Both experimental and analytical 
investigations are reported. Three different analytical models have 
been formulated and are described in detail. An experimental 
investigation has been conducted on multiple grooved two-inch diameter 
viscoseals over a wide range of gas densities and shaft speeds up to 
30,000 rpm. Comparisons 'are presented between actual viscoseal 
performance and the theoretical predictions for both sealing coefficient 
and net leakage parameters as functions of the degree of gas rarefication. 
Recommendations are presented for the use of the analytical models. 
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NOMENCLATURE 

Groove , aspect ratio, b'/h 
Hypothetical groove aspect ratio 
Groove cross sectional area 
Rotor induced flow coefficient 
Annulus flow coefficient 
Dimensionless annulus flow coefficient 
Continuum limit of A^ 

Manometer cross section area 
Groove pressure flow -coefficient 
Rotor induced flow coefficient 
Groove axial width 
Groove normal width, b cos a. 

Groove flow coefficient 
Radial seal clearance 
Mean radial clearance , 

Shaft diameter 

Slip coefficient constant of proportionality 
Groove depth 

Manometer total deflection 
Manometer -deflections 
Boltzmann constant 
Seal radius ratio s 1 - 2c/d 
Flow passage length coordinate 
Tube length coordinate 
Groove length coordinate 
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Land length coordinate 
Total groove , length on seal 
Seal axial length ■ 

Sealant mass per molecule 

Sealant mass flow rate 

Sealant molecular weight 

Molecular density, N/V 

Number, of thread lands in seal section 

Number of molecules i 

Sealant 'molecular flow rate" 

Seal speed, rpm 

Knudsen number, A/c 

Groove Knudsen number, A/b f 

Knudsen number based on groove depth, A/h 

Number of groove starts 

/ 

Total turns of spiral on seal 

Absolute pressure 

Average seal pressure, (P^ + 

Atmospheric pressure , 

Upstream seal pressure 

Downstream seal pressure 

Differential pressure across seal, P^, - P^ 

Sealant volumetric flow rate 

Volume flow rate at unit pressure, V'P 

Volume flow rate at unit pressure, in groove 

Volume flow rate at unit pressure over land 
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Net volume flow rate at unit pressure 

Rotor induced volume flow ,rate at unit pressure 

Inner radius of annulus .(shaft) 

Outer radius of annulus (housing) 

Seal pressure ratio, or ^ 2^1 

Distance ratio, r/b f 

Distance ratio, rr/b T 

Universal gas constant 

Tangent of helix angle a 

Flow measurement time increment 

Absolute temperature 

Velocity components in £, n » z coordinates 
Average rotor induced velocity 
Surface velocity 

Dimensionless, groove flow velocity 
Effective rotor velocity 
Groove flow velocity 
Rotor induced flow velocity 

Rotor induced flow velocity in hypothetical groove 
Volume 

Mean molecular speed 

Volume flow 'rate 

Total groove volume flow rate 

Volume flow rate in a single groove 

Continuum limit. of V' 

g 

Land leakage volume flow rate 
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Total volume flow rate in seal 

Land axial width 

Land normal width, w cos a 

Manometer fluid specific weight 

-1 n s (b ’ + - w *) 

Seal helix angle, sin [ ^ ] 


Eigenvalue 

Seal geometric parameter, (h + c)/c 
Tube flow coefficient 
Seal geometric parameter, b/(a + b) 

Tube flow coefficient 
Coefficient of slip 
Mean free -path 
Sealing coefficient 
Absolute viscosity 
Pressure unit, micron of mercury 
Density 

Ratio of .molecular collisions 
Elemental volume 
Rarefied viscoseal parameter 
Solid angle * 

Angular velocity 
Special Subscripts 
r Refers to land 

g Refers to groove 

£r Refers to £ direction in region of land 
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T|r Rerers to n direction in region of land 

£g Refers to £ direction in region of groove 

rig Refers to r) direction in region of groove 

y Refers to axial component 
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1. INTRODUCTION 

An investigation of rarefied, internal gas dynamics with emphasis 
on shaft sealing applications was initiated on May 1, 1966 at the 
University of Tennessee in the Department of Mechanical and Aerospace 
Engineering, The investigation was conducted for the National .Aeronautics 
and Space Administration under Research Grant 43-001-023, This report is 
the final report of this investigation with major emphasis on the 
rarefied-gas viscoseal. 

The initial efforts were directed at a number. of basic problems. 

The results of these investigations have previously been documented 
and only the references will be repeated here. These investigations 
concerned rarefied flow through short tubes (1,2) , annuli (3,4), long 
square tubes (5), and other basic geometries (6-8), From these 
investigations a more complete understanding of rarefied flow has 
provided the basis for analysis of more complex applications. 

The major effort, of this study has been directed to developing 
analytical models which can be used to predict the performance of 
viscoseals operating with a rarefied gas as the, sealant. Since no 
experimental data were available which were adequate for establishing 
the validity of analytical models a significant portion of this effort 
has been devoted to the development of such data. This report details 
the analytical models which were developed and also the experimental 
procedures and results. Conclusions are presented based on the 
comparison of results from analytical models, with experimental results." 

2. ANALYTICAL MODELS 

Three different analytical models have been developed in this 
investigation and will be presented in detail. The first is based on 
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Hodgson’s (9) work and is called the modified Hodgson Model. The next 
model is based on the work of Boon and Tal (10) and is called the Slip 
Modified Boon and Tal Model. The third model is based on a super- 
position of the groove and annulus flow components and is called the 
Annulus-Groove Model. 

2.1 Modified Hodgson Model 

In Hodgson's (9) investigation only seals with a single thread 
start were considered. The development presented follows closely that 
of Hodgson but the analysis has. been generalized to include multiple 
thread start seals and also a technique is presented for solving the 
specific molecular flow rate equation for a sealed condition. 

Hodgson chose to analyze the particular configuration where a 
smooth shaft rotates within a grooved housing. It will .be shown later 
that -the analysis for the grooved shaft and smooth housing is Identical 
to this configuration. This being the case, HodgsonJs model is not as 
restricted in this respect as it might ■ appear. 

Hodgson considers the flow dn the viscoseal to be composed of 
three basic components: (1) the pressure induced flow along the groove, 
(2) the pressure induced flow over the lands, and (3) the rotor induced 
flow in the groove. This treatment of the flow is quite common and is 
exactly th,e way King (11) chose 'to divide the flow for analytical 
considerations . Figure 1 shows a representation of the flow components , 
Q^» Q^> and Qg represent the pressure induced land flow, the rotor 
induced groove flow and the pressure induced groove flow, respectively. 
Hodgson also restricts his' considerations to a -seal with a single thread 


start. 
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In his development of the two pressure induced flow components , 
he takes the semi-empirical approach of Knudsen (12) in describing the 
flow throughout the entire regime from continuum to free-molecule flow. 
In his long tube work Knudsen found that his experimental data, could 
be described by an equation of the form 


Q' 


[y'P + 6 


, i + vp i 

1 + £p J 


dP 

di* 


( 1 ) 


dP 

In the continuum limit Knudsen reduced Equation (1) to Q’ = y'P 
and by equating this to the known continuum. solution, he was able to 
determine y’. In a similar manner as the free-molecule limit was 
approached he reduced Equation (1) to <*’ = **§ . and by equating this 
to the known free-molecule solution .he determined 3’. By considering 
the slip flow regime, Knudsen was able to determine the ratio V/£' 
by noting that Equation (1) becomes Q' = (Y ? P-+ and equating this 

to the known continuum with slip solution. Next Knudsen determined 
the difference £’-V from a consideration -of nearly free molecule flow 
and thus was able to determine both and V„ He then applied experi- 
mentally determined corrections to £' and V so that the experimentally 
observed minimum in the Q*/AP versus P curve would be correctly predicted. 
Knudsen .applied the analysis above to a long circular tube. In his 
analysis Hodgson applies the identical procedure to his treatment of 
the pressure induced flow in the groove (a long rectangular duct) and 
the pressure induced flow over the lands which he takes to be a narrow 
slit. 


2,11 Pressure Induced Flow in the Seal Groove . Hodgson assumes 
that the pressure varies continuously along the axis of the seal groove 
from a value of P 2 at the high pressure end to P^ at the other 
end. In order to determine the pressure gradient 
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gradient along the ’groove axis, it is necessary to relate the length 

of the groove, $> , to the seal length, L. A development of a viscoseal 
S 

is shown in, Figure 2. Hodgson only considered seals with a single 

thread. The development that follows is generalized to any number of 

thread starts, n . 

s 

In Figure 2 line: AB is drawn perpendicular to the grooves. The 
number of turns of spiral that AB crosses is equal to the total number, 
N fc , of complete turns on the seal which can be expressed as 


„ _ h cos a 

t w* + b’ 

The length of groove per turn of spiral is TTD/cos.ot, The total groove 
length on the seal is 


Z 


T 


TO L cos d _ L7TD 
cos a'w’+b' w r + b ' 


The length of each groove is then 

n _ &T _ LTTD 

g N N (w ' + b ' ) 
s s 

It follows that 


dP dP 

dH ~ LTO , 

8 L N (w t + b r ) J 


N (w* + b 1 ) 
s 

TO 


dP 
dL * 


Hodgson -indicates that the groove pressure gradient in a single threaded 
seal is 


dP = w’ + b T dP 

d£ TO dL 

g 

from which it follows that the flow in a groove of a multi-threaded 
seal is N times the flow in the groove of a single threaded seal of 
the same groove width, land width and diameter. Since there are N s of 
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2 

these grooves, the total groove flow in the multi-threaded seal is N 

s 

time the flow in ,a comparable single threaded seal. The total groove 
flow is then 



= N 2 [(BP + C 
s 


1 * C 1 P dP, 

1 + C 2 P } dL J * 


( 2 ) 


where the bracketed term is the single threaded groove flow developed 
by Hodgson. The constants B, C, C^, and C 2 depend on the geometry of 
the seal and the properties of the sealant and are given in Appendix A. 
Appendix A also contains a discussion of the flow models used by 
Hodgson to obtain these constants. 

2.12 Pressure Induced Flow Over the Lands . Based on the assumption 

that the pressure varies continuously along the helical groove, Hodgson 

shows that the effective pressure gradient for the land flow is 

dP w* + b' dP 

dx-^ w dL 

This pressure gradient applies .equally well to both the single and multi- 
threaded seal. The land leakage flow for a seal of any number of threads 
is then 


1 + C 3 P dP 

t D ' p+I TTc?p]§ 

4 


(3) 


where the constants D 1 , E, and also depend on the seal geometry 
and the properties of the sealant and are given, in Appendix A along 
with a discussion ,of their origin. 

2.13 Rotor Induced Flow. Hodgson takes a simplified approach to 


the prediction of the rotor induced flow. In Figure 3 a groove cross 
section is shown with .the rotor .moving over the top of the groove. 
Hodgson develops the rotor induced flow on a molecular basis, but as he 
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points out a continuum, approach yields the same result. It is assumed 
that the rotor induced flow is the same as the flow obtained in a long 
rectangular duct of width b ' and height h in which the upper-wall moves 
with velocity U cos a (the component of the circumferential velocity of 
the rotor along the axis of the groove). Rather than solve the describing 
differential equation for parallel flow, Hodgson chooses to compute 
the volume flow .based on an area-weighted average velocity. This average 
velocity is 

~ _ (U cos cpb T + (2h + b 1 ) (0) Ub * 

u 2 (b ' + h) ~ 2(b f + h) cos a ° 


The total rotor induced flox* is thus 


Qr 


2 

fe ) c ° 8 al P - H S A ' P 


(4) 


where is generalized to a multi-threaded seal. 

2.14 Total Seal Flow . The flow rates given by Equations ’ (2 , 3, 
and 4) are superimposed to give the total flow in the seal which can 
be expressed as 


2 ^ 4" , 1 + C P . 

% - - N S [BP + c TTcT 3 ^ - [D-P + E —fa f - N S A’ P < 5 > 

2 4 


where flow ,in the direction of decreasing pressure is considered positive. 

2.15 Solution of the Modified— Hodgson Equation . Unlike Hodgson’s 
basic equation. Equation (5) is applicable to a seal with any number of 
threads. There are three solutions to Equation (5) which are of 
particular interest. These three .cases are:. (1) the flow rate through' 
a non-rotating seal, (2) the flow rate through a rotating seal^ and 
(3) the pressure , difference across a rotating seal when the net flow 
is zero. 
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Flow Rate in a Static Seal 

The volume flow rate at unit pressure, Q’ = V ? P, is related to 

the molecular flow rate, n, by 

*01 

n ‘ kT * 

For the case of a non-rotating seal (Q^ = 0) , Equation (5) can be 
integrated to obtain the specific molecular flow rate 


4 0 C 

Ip-SL twjB + D')P + M s 2 C c7 +E 07 

2 4 


0 C. - C- r + 1 r (1 + 2C„P) + 1 

+ N_ C — (■ / ■ _ ■ y) In [-z 

p r + 2C 0 P + 1 


20 2 P 


] 


+ E- 


C, - C. r + 1 r (1 + 2C.P) + 1 

4r— 1 r> m t-2 


2- 
2C7P 
4 


3> 


( 6 ) 


r + 2C.P + 1 
P 4 


Flow Rate in a Rotating Seal 

Assuming a constant pressure gradient along the groove, the total 


rotor -induced flow on a, molecule basis is 

rh 


n 

ap 


N A’P 
s 


N A’ r 4- 1 

S [— E. 


r - 1 
P 


(7) 


kTAP 2kT 

The net flow through a rotating seal is found by combining the rotor 

induced flow of. Equation (7) with the flow in the static seal. Equation 

(6). The net, specific .molecular leakage is 

C 


i?" hl { < 8 s 2b + d ') p + ^ c 5- +e c7 

2 - 4 


r + 1 




C_ - C 1 r A (l + 2C.P) + 1 

2 1 ' ' 1L 2 .> 


2 

C 2 P 


r + 2C_P + 1 
P 2 


C 4 " C 3 , 

+ E - - • • In (■ 

C 4 J 


r (1 + 2C.P) + 1 
4 

r + 2C.P-KL 
P 4 


■) - N g A'L] }. 


( 8 ) 
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Pressure Difference at Zero Net Leakage 

Since a rarefied viscoseal, in the ideal case of a true space 
environment, will normally operate with P^ = 0, it would be impossible 
to maintain a zero net flow. However, the condition of zero net, flow 
is of interest as far as experimentation and comparison to continuum 
performance are concerned, and in non-space applications where P^ 0. 

The maximum pressure difference under which a seal can maintain 
a zero net flow is found by solving the modified-Hodgson model equation 
subject to the condition that = 0. Integrating Equation (5) subject 
to this -condition one obtains 


2 *"”1 ^3 

2P + AP -,N^C+E 2 + C 2 (2P + AP) N s C<: cJ _1) 2 + C 4 (2P + AP) E ^ " 1 ' ) 


] g 

2P - AP 


2 + C 2 (2P - AP) 


2 + C 4 (2P - AP) 


= .exp [fo A’L - (N B +D l )AP], 

{ 3 S 


(9) 


In general Equation ’(9) cannot be solved explicitly for AP. For the 
special case of continuum flow 2P > > AP,the continuum AP can be 
obtained from Equation (9) as 


AP = 


N A’L 
s 


n B + D 
s 


( 10 ) 


At the other extreme, free-molecule flow. Equation (5) can, be reduced 
to its free-molecule limit. 


On = ' CN s c + E >f - V p . 


(10a) 


Since in the free-molecule limit P is very small, the last term in 
Equation (10a) could presumably be > very small; but since efficient seal 
performance requires that the rotor induced flow be of the same order 
of magnitude . as the pressure induced flow, this term is retained. 
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Solving Equation (10a) subject to the , condition that =,0, the free- 
molecule pressure ratio becomes 

r = exp[N A'L/(N 2 C + E)J. (11) 

F S S 

The AP across the seal can be expressed as 
r - 1 

• < 12 > 

P 

Combining Equations , (11) and (12) , the AP in the free-molecule regime 
is 

exp [N A'L/(N 2 C + E)] - 1 

AP = 2P - ^ . (13) 

exp[N A*L/(N C + E) ] + 1 
s s 

Since no explicit solution to Equation (9) can in general be 
obtained, some approximate solution technique must be employed. Since 
Newton’s method of approximating roots is generally a -rapidly converging 
iterative method, it is employed in solving Equation (9), Newton’s 
method requires that the given relationship be . differentiable , 

Equation (9) can certainly -be .differentiated. Another very important 
requirement is the ability to, make a close initial approximation to the 
solution. This requirement is particularly important with a complicated 
relationship such as i Equation (9). 

In order to make a close approximation , to the roots of Equation -(9) , 
one needs -to know as. much as possible before hand about the . character 
of the solution. At this point two characteristics are known? (1) the 
continuum limit, Equation (10), and (2) the free-molecule limit; 

Equation (13) . The solution .of Equation (9) will simply define the 
behavior over the entire range of average pressures and approach the 
continuum and free-molecule limits at the two extremes » In order to 
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make a close approximation to the value of AP which satisfies Equation 
(9) for a given, P, a high average pressure, P^, is chosen initially such 
that the continuum solution from Equation (10) is ?a good approximation 
to the root of Equation (9). Using this initial approximation for AP, 
Newton's method is employed to solve Equation (9) for the AP at P^« 

An incremental decrease, AP, in the average pressure is then taken 
and the AP obtained from the .previous iteration at P^ is used as the 
initial estimate of the solution for AP at and the, iteration process 
is repeated to obtain AP. The entire process is repeated ,to obtain 
AP at P 3 , P^, etc. 

Essentially the above process could be continued until the entire • 
spectrum of pressures had been traversed, but the process encounters 
•difficulties' near the intersection .of the continuum and free-molecule 
asymptotes. In the near continuum region, the AP versus P solution 
of Equation (9) is fairly flat, thus -making the solution at ^ a 

good approximation to the root at P^. In the intersection region, 
however, the rate of change is -so large that the- method used above 
for the initial estimation 'of the solution is not sufficiently accurate. 
Two simple modifications to the above method can help to, ensure a 
close approximation of the root at -P^. An obvious modification would 
be to reduce, the step size, AP, Since all previous points P^ through 
P^ _ ^ have been determined, these points can be used to extrapolate 
to an initial approximation ,at P . Essentially an extrapolation of 
order n - 2 could be made , but experience has shown that a linear 
extrapolation combined with successive reductions in AP is sufficient 
to ensure convergence at P . 
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After passing through the critical region in. the vicinity of the 
intersection of the asymptotes , the curve essentially assumes the 
straight line predicted by Equation (13). The solution is thus complete. 

The Fortran programs for the solution of Equation (5) for the 
three cases outlined above are presented in Reference 13. The flow rate 
in a static seal is simply a special case (Q^ = 0) of the rotating 
solution. Consequently only one program is needed for the 'flow rate 
solutions , 

The modified-Hodgson model provides a convenient means of predicting 
the performance of a given viscoseal. The limiting conditions of 
Equations (10) and (13) can conveniently be, used to predict optimum 
seal geometries in .the purely continuum and purely free-molecule 
regimes, respectively. In the transition regime, however, Equation (9) 
would be very difficult to use in an optimization study because of the 
time consuming solution method. 

2.2 Slip Modified Boon and Tal Model 

In 1959 Boon and Tal published (10) a significant analysis of the 
viscoseal in the laminar continuum, regime for a constant density fluid. 
The viscoseal geometry was ’approximated by two flat plates, one of which 
was grooved, moving parallel to each other. The flow .field was developed 
from a superposition -of Couette flows along and across, the grooves. 

The resulting velocity distributions were integrated to obtain the 
volumetric flow rate, and the pressure generation , for zero net leakage 
conditions was developed. Subsequent investigations by Stair (14) 
using liquids as the sealant showed very good agreement with the analysis 
of Boon and Tal. The data of Hodgson and Milligan (15) using air as the 
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sealant' also substantiated this .analysis for laminar continuum operation , 
with gases. 

The Slip Modified Boon and Tal Model has been developed by formulating 
"rarefied" corrections to the laminar continuum equations , .developed by 
Boon and Tal. This approach is similar to that used in previous 
investigations (3, 4, 16, 17) of internal rarefied flow that have 
indicated that a .single model for gas flow through tubes 'and annuli, 
applicable to the flow regime extending from continuum to free molecular 
flow can be derived by the inclusion ( of rarefied effects on the continuum 
model. 

Consider a screw formed on , a. shaft located concentrically within-a 
cylindrical .housing with a radial clearance c„, The annular space is 1 
filled with a gas and the shaft is moving relative to the .housing with 
an angular velocity, Figure 4 shows a developed view of the 
viscoseal geometry. The (x, y) axes are along and -normal to the 
direction of relative motion and the (£, T]) axes are parallel and 
normal to the grooves. The (x, y) and .(£, -r|) coordinates -systems are 
related by: 

£ = x cos a + y sin a 

n = y cos a - x sin a. 

Previous investigators (10, 18) have reduced the describing partial 
differential momentum equations to the Reynolds lubrication equations. 

These -analyses were based on the .flat plate model of Figure 2 and 
assumed steady, isothermal, two-dimensional, laminar flow. ‘Further, 
these. analyses .assumed inertia forces to be negligible in comparison to 
viscous . forces and neglected body forces and end effects. The 
describing mathematical model is taken as: 
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d 2 u = 1 3P 
dz 2 " ^ 8 « 


(14) 


d 2 v = _1 3P 
dz 2 “ ^ ^ * 

Integration of (14) and (15) gives: 


(15) 


u 


1_ 3P 2 
2U 9? Z 


+ Cz + C 2 


(16) 


v 


3_ 3P 

2y 9n 


2 

z 


+ C„z + C. . 
3 4 


(17) 


The integration constants are determined by the boundary conditions. To 
account for the non-continuum effects, slip boundary conditions are 
introduced as derived by Kennard (19), For a monatomic ,gas , flowing within 
two parallel boundaries , the slip boundary condition can be written 
in the following, general form; 

u | = u - - + GX (■!“) | n 

gas'z = 0 wall 3z ‘z = 0 


+ !^i>U-o + °a 2 > 


(18) 


where u is the velocity in .the tangential x direction, z is th.e 'direction 
normal to the motion, and G is a proportionality constant. If the 
temperatures of the surfaces are assumed constant and equal to the gas 
temperature then the 'Slip boundary condition reduces in ..general form to 


u = u . - + GX (-lr~) = ,0, 

gas wall da z 


(19) 


For the particular geometry under study, the slip boundary conditions 


are: 


Along the lands: 


du 


u t „ = U cos a + GX — 5 — „= u- 

r z = 0 dz z = 0 1 
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dii 

U r ^ z = h dz ^ z = h “ U 2 

r r 


Along the groove: 


dll 

u I n = U cos a + G\ - 5 — I 
g'z = 0 dz'z-o 

da 

u gl z -h ”- GX d^l z =h E V 

g g 


= U, 


Across the lands: 


r 1 z - 0 


= 0 = - U sin a + GX ~-\ 
dv 


dv_ 

dz ‘I z, = 0 = V 1 


v rU = h G \lz ' z = h v 2 “ 

r r 


Across the groove : 


dv 


v | _ n = - U sin a + GX [ _ = v„ 

g'z =0 dz 1 z = 0 3 

dv 

v I , = - GX -?-£■ | = v. . 

g'z=-h dz 1 z = h 4 

g g 

Using the above boundary conditions the following velocity 
are determined: 


Along the lands: 
1 3P , 2 


u r 2 y 3£ 


u - u. 

(a‘ " h r z) + (■ - h --) z + u r 


Along the grooves: 
1 3P , 2 


U g “ 2 y 3C 


u, ~ u_ 

(a'"' - h z) + (-— -r — — ) z + u,, 
g ti g 3 


Across the lands: 

1 3P , 2 , N , , v 2 “ V U 
v r ‘ 2 y 3? <Z " h r 2) + ( — )z * v V 


Across the grooves: 

1 3P , 2 . /4 - V 3, , 

V 8 = 2?'ST (a - V 5 ( “ )z+v 3' 

o 


components 


(20) 


( 21 ) 


( 22 ) 


( 23 ) 
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The slip velocities at the walls can now be determined, and the following 
velocity distribution, equations obtained: 


1 8P , 2 , , , ,-U cos a . z 

u r “ 2]J 85 (z " h r z ^ + ( 1 + 2GA/h r ) h r 

1 + GX/h r 1_ 3P . 2 GA 

+ U cos a ( 2GX/h ) ~ 2]1 3^ h r h 

r r 

1 8P , 2 . v - ,-U cos os z 

u g 2y 8t (z " h g z) + ( 1 + 2GA/h ^ h 

g g 

1 + GA/h , 3 p ? p)i 

+ D 008 “ c rr^7h > - 2U r 

g g 

1 8P,2 , \ sin a z 

v r 2y ‘ 8u (z h r z) + ( i + IGX/hJ h r 

it • 1 + GX/ \ 8p 2 GA 

U Sin a( l + 2GA/h ^ 2U ^ h rh 

r r 


v — 1— — (z? - h z) + s - - a — ) — 

g 2y 8n U g J YL + 2GA/h ; h 

g g 

• 1 + GX/h g __ 1_ 3p 2 GA 

U sin a( 1 + 2G X/h ^ 2li 3n h g h * 


Noting that the axial velocity components are: 

u^ = u.sin a (28) 

Vy = v cos a (29) 

the axial flow rate components and Q^g be , determined. 

The width of the flow path for the 5 land flow component is (1 - y)ttD 
and the path width for the 5 groove flow is yTrD. The. ratio of the groove 
width to the groove plus land width is defined as y„ The axial 
component of the 5 coordinate land flow is: 


h h 

Qey = (1 ~ Y)T® / r U dz = (1 - Y)TTD / r U sin a dz. 


Substituting Equation (24) into Equation (30) and integrating gives 
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Q 5r = (1 - Y)7TD sin at— (||) r h 3 (- - -■) + U 'Cos Ot~ ]. (31) 

In a similar manner: 

hr B sin 4f> g h g r - - i> + D cos ° ^ < 32 > 

g 

V - a - Y>™ COS «[ij;(||) r h 3 (--S. - £) - U sin OjE] (33) 

ll 

%, - ym cos afij (f) g h 3 <- - ±) - D sin a fl . ■ (34) 

s 

The pressure gradients in Equations (31) through (34) may be replaced 
by the more convenient axial gradients by noting that 


3P 3P . 
r? = -k— sin a 
3y 


and 


(1 


-™f> r 


c iE.) 


„ is. 

g 3 y 


cos a 


(35) 


(36) 


From the continuity of mass, the flow across the groove in the r| 
direction must equal the flow across 'the land in this direction „■ Thus 
it can be shorn that 


p l / 9F- . j3 , GA 
L 2 y t '3n ; r r h r 

,1 )3F. 3 . GX , 

2y Wg h g ( ~ h " 

© 


- t) 




h 

- U sin a ~ 3 = 

h 

- U sin a ]. ■ 


The total flow is given by 


(37) 


Q = Q ?r + Q £g + V + V C38) 

Using Equations (31) through (38) , the expression for the axial flow 
is : 


1 9P 
12\i 9y 


c 3 7rD 

11 + t 2 ) 


[t 2 a - Y)d + ~) + a +;•—■) (3 3 )(Yt 2 +|)i 


Q = - 
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4 * 

UCTTD t(l - Y)(P ~ 1) r P U Be J 


~ i + t 2 % a + l£i) ' 


where 


B = 

Y “ 


h + c 
c 

b 

w + b 


1 + WX 

f - (1 - Y)B 3 ( fgjj— > + Y 

+ ’c 

t = tan a 

The Knudsen number for the viscoseal is defined as 


H =A 
V c » 


( 39 ) 


(40) 


where seal radial clearance* c, has been selected as the characteristic 
system dimension because of the manner of its appearance in Equation 
(39) in relation to the mean free path A. 

Equation (40) can be rewritten ‘in terms of molecular flow rate 
rather than volumetric flow rate as 


* P 3P C'TTO r .t.2 /t , 6GA^ , 6GX x / d 3 n 

» " - a? 77772 , [t a - v)(i + — ) + a + er><s > 

\1> T t J ry fipX 

(Yt 2 + 1)] + sp- f a + vr ) - lh (41) 

* 2 l + t 2 Hi+^) 


Further consideration now needs to be given to the slip coefficient 
constant of proportionality "G". Various investigators have developed 
expressions for this proportionality constant ranging from 2/3 (20) 
to (2 - f) /f (19) where f is the fraction of their tangential momentum 
which molecules give up upon striking a solid boundary. Published 
values of f by Kennard (19) range from 0.79 for air flowing over fresh 
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shellac to 1.00 for air or C0 2 over machined brass. A value .of f 
equal to one. is taken to be realistic for the flow of gases over 
machined surfaces. Thus the value of the slip coefficient G was 
taken to be unity. 

Knudsen n umb er as defined in Equation (40) indicates the degree 
of rarefication since -it becomes large as the’ mean free path, X, 


becomes large. When the Knudsen value becomes .small and approaches zero, 
corresponding to continuum conditions, then Equation (30) reduces to 
the continuum, no slip, solution shown by Stair (18). 

For isothermal flow the mean free path, X, is a function of the 
pressure and gas properties., The Chapman relationship (19) can be 
used to express this dependence. 


X = 


V (2& 1/2 

-TO , 

2p 


(42) 


This relationship -will be ,used to express the right hand side of 
Equation (4l) in terms of pressure. 

2.21 Variable Mean Free Path Solution . In order to solve 
Equation (41) for the flow rate in terms of the seal length and the Ap • 
across the seal a numerical scheme must be employed. Before integrating 
the equation is rearranged. * Let P q be the .pressure at which 
X = C/6G., Then 

uaA 112 

6GAP 6G UU7r ’m ; 


Now Equation , (41) may be written in terms of pressure 


P 2 P 1 c 3 7TD 


n = - 


dP' 


2 dy [t z (i - Y)(1 + fr) + a + gjrr)3 (yt + ^)] 


12ykT(l + t) 

+ ueP o P'7rDt(l - y)(3 - 1) ' 8 3 (1 +■—,) 

' “ o t i “ 


2kT (1 + O 


^(1 + | r ) 


( 43 ) 
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where 


P* + P/P 


and 


* a + ftlr) 

$ = (1 - Y)3 r— + Y 

(1 +~r) 

Non-r-dimensionalizing the length coordinate, y, by using the seal 
length, L, i.e., y' =-y/L and dy = Ldy* we can write 
n = f(P') g(P') 


where 


f(P') = - 


2 3 

P C 7TD 
o 


12ykT (1 + t ,)L , 


[t z (l - Y) (P* + 1) 


(44) 


+ (P* + ^)g 3 (Yt 2 + |)] 

cUP 7TOt(l - y)(3 - 1) o CP’ + 1/3) 
g (P T) [3 3 - P'3 

2kT (l + O ifj(l + 1/p') 

Simplifying Equation ,(44) by defining a new flow rate parameter n.we 
obtain 

7" jpi 

n = f dp" + A f g (45) 

where 


A = 6 ^ UL 

f 2 

c P 

o 


f - f(P' , 3 ,Y»t) 


1 [t 2 (l - Y)(P’ + 1) + (P l + j) 


(1 + O 


3 3 (Yt 2 + |)3 
g « g(P' ,3,Y.t) 


e t(i - y)( 3 - l> 

(1 + t 2 ) 



< pl + 

m + jr) 


p'3 
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Separating variables in Equation (45) and integrating over the 
seal length gives 

1 = / 2 _ dP* (46) 

n - A f g(P') 

For a given physical condition,, i.e., A^, the problem 

reduces to solving Equation (46) for n. Newton's method was used by 
letting 

r P 2 3 P,) dP' = 0 (47) 

' 1 " J p , ft -Ag(P') 

Newton's method predicts successive approximations to the solution of 
Equation (45) through the following relationship: 


n , = n - F/F r In 

n + 1 n 1 

where 

r - ^ - ! 2 -|£h_ , dP- 

d A PJ (n - A f g (P ' ) ) Z 

Since the most useful situation for the viscoseal as a sealing device 
would be to have the net flow rate equal to zerp a similar technique 
was used to solve Equation (41)i for the AP which results ,for n ,== 0. 

A widely used parameter in sealing work is defined for no net flow 
in the form of sealing coefficient, A. 


_ 6tiDL 
? 

APc 


(48) 


Results for the variable mean free, path will be presented for both net 
leakage and for no net leakage (sealing coefficient;) values. 
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2.22 Mean-Value Solution . If the mean value of seal pressure is 
used in Equation (42) to calculate a mean value of mean free path, 
then Equation (39) can be approximated by replacing the local mean 
free path by an average value. This permits the simple integration of 
Equation (39) to give 


Q = " lk ° ™ 2 £ 1<V + g gg) 2- - f 2 <V 

12li QL (1 + t 2 ) 1 K 2(1 +t 2 ) 2 K 


where 


OvjriN O 1 

f x CN K ) - t Z (l - Y)(l + 6GN r ) + (1 + -™)(3 J )(Tt + 


f 2 (N K ) = t(l.- y)(B 


- i) L 


$ 3 (i + 6G \) 

[ g 

i^(l + 6GN k ) 


For no net flow the sealing coefficient can be written as. 


A = 


6]iUL 

APc 2 


f l C V 

f 2 (N K ) 


(50) 


Results will be presented using' this mean-value type solution. 

2.23 Particle Slip Correction . In the previous section where the 
slip-modified Reynolds solution was derived, slip boundary conditions 
were applied to account for the decreased intermole cular momentum 
transport at the walls. The analysis attributed a slip velocity -relative 
to the walls .for all the molecules adjacent to these surfaces. For 
the case of stationary walls, this implies that every molecule, on the 
average, will possess additional flow velocity due to the slip flow 
contribution. The preponderance of evidence presented by Weber (16) , 
Kennard (19), Present (20), and Fryer (21> indicate that molecules 
whose last collision was at the ; wall can have no slip velocity since such 




- 22 - 


molecules are diffusely reflected from the walls. Thus only those 
molecules coming from collisions with other molecules can possess a 
slip velocity. A correction will now be developed to account for those 
molecules which do not experience slip boundary conditions due to their 
collisions with the walls. The correction .will be established by 
determining the ratio of .the molecule-to-molecule collisions to the 
total number of collisions, or the sum of the molecule-to-wall plus the 
molecule-to-molecule collisions. This ratio, a, represents the 
fraction of molecules present which experience the slip boundary condition. 

The number of molecules .striking the wall per unit time and unit 

n y 

area is (20) — where V is the mean molecular speed. If >one considers 
the unwrapped viscoseal geometry of Figure' 4, the seal surface area is 
approximated by 


27TDL + 2h ra „ (— Jj-r-) 

cos a s w + b 

where the right hand term, is the surface area contributed by the groove 
side walls. Thus the number of molecule-to-wall collisions in unit 
time is 

— ^ nDL[l + > — . 

2 cos a(w + b) 


The total number of molecule-to-molecule collisions per unit time 

xiV 

and unit volume is (20) . The volume of the viscoseal is the 

volume of the .annular space plus the volume of the grooves . Thus the 
total number, of molecule-to-molecule collisions per unit time is 
nV r „ T ,2 2 n , 7lDL*h’b , 

— [,,L(r 2 - r i> + 15TT) 1 - 



(r 2 + rp(r 2 - rp , 


Since 
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and, due to the small clearance of the normal viscoseal geometry, 
'u 

r^ + r^ h* D and - r^ = c, 
then the expression simplifies to 


hV -n-Tvr r , h • b, 

T~* TOL [c H T-r" 

A w + b 


Hence, the proportion of the molecules which experience slip .is 


0 = 


1/N r [1 + y(B - 1)] 

1/n k [1 + Y(3 - 1)]+ -|[1 + 


cos a(w + b) 


with 

c 


X = 


b 

w + b 



and 



In the solution of the .slip-modified Reynolds equation, it will be 
recalled that boundary conditions were not permitted on the groove side 
trails .. Thus it is more in keeping with the manner of this analysis to 
omit the molecule-to-wall collisions .which occur on these surfaces 
and contribute the h/cos a(a + b) term in the .denominator of Equation 
(51). With this restriction, the equation becomes 


1 + N k / 2[1 + (3 - l)y] M 

In the analysis of slip flow through long tubes by Weber (16) , 

Fryer (21) , and Milligan (17) , the resulting flow equation could be 
manipulated into two separate terms. One of these terms is the . continuum 
Poiseuille flow term while the second term is the slip flow contribution 
which results from .application of the slip boundary conditions. 

Equivalent, type terms were obtained by Milligan, Cowling, and 
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Wilkerson (4) in the analysis of rarefied flow in a long- concentric 
annulus. For both of these geometries » the correction for the particles 
which do not experience slip was made by multiplying the ,slip flow 
contribution by C, since it is this uncorrected term which implies that- 
all the molecules possess the additional slip velocity. Thus when.cf is 
applied in this manner the resulting slip contribution term has been, 
corrected for the particles which experience wall collisions and thus 
possess no slip velocity. When this correction is made to Equation ,(49) 
the following equation is obtained: 


Q 


1 AP 
12V L 


c 3 ttd 
(1 + t 2 ) 


W + 


UcTfD 

2(1 + t 2 ) 


vv* 


(53) 


where 


f 3 (N K > = crf^NjP + c 5 (l - <0, 

and 


f 4 (N R ) = crf 2 (N K ) + c 6 (l - cr) , 

r = Yt 2 (l -Y) ( 3 3 - l ) 2 + e 3 (l + t 2 0 
“ 3 

fT(l r Y) + Y 

c = Yt(l ~ Y)(g 3 - 1)(3 - 1> 

6 g 3 (l - Y) + Y 

In a similar manner to that 1 of 'Equation (50 ), the sealing 
coefficient for the mean-value .slip-modified Reynolds solution with 
corrections for particles which do not experience slip is 


A 


6yUL . 
A Pc 2 


f 4 (N K ) * 


' (54) 
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Often it is more desirable to express results in .terms of molecular 
flow rates rather than volumetric flow rates. Equation (49) for the 
Mean-Value Slip Modified Boon and Tal solution with no correction due 
to particles which do not experience slip' is 


n ' 
AP 


, it 1/2 , C TTD 
'■2mKT ; „ . J2 


\ (1 + t )L 


) f x (N K ) 


r + 1 

+ (rf— r) 


UcTTD 


p 4KT (1 + t ) 


W> 


(55) 


where the relationship that 
_ r - 1 

AP = 2P ' ( - r - P - T y ) , (56) 

P 

has been written in. terms of the pressure ratio s r , across the seal in 
the right hand' tern. . Similarly* the specific leakage flow rate with 
the correction for the particles that do not experience slip is 


n . TT . 1/2 c 2 TTD 

AP - - W ( + c 2 )l 

lx 

r + 1 

+ ( _E 

r p " 1 ‘ 4KT(1 + t 2 ) 


-) f 3 (N K ) 

W" 


(57) 


It should be noted that for a fixed .seal geometry .and a. given gas 

at a specified temperature, the specific molecular ,flow .rates of 

Equations (55) and (57) are 'functions only of Knudsen number, the 'speed' 

U, and the viscoseal pressure .ratio, r . As -the pressure ratio becomes 

high, then (r + l)/(r - 1) approaches unity and the specific flow 

P P 

rate is a function of the speed and the -Knudsen number only. 

This same technique can be used to evaluate the particle slip 
correction for the variable mean-free path solution by correcting 
Equation (43) before using a numerical scheme to solve the , equation. 
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2.24 Self-Diffusion . The terra diffusion, as used here, refers to 
the molecular transfer which occurs due to a concentration .gradient. 

If the gas is a pure unmixed gas, then the diffusion is one of self- 
diffusion. In their work concerning low density flow of gases in a 
capillary, Pollard and Present (22) suggested that at .low pressures the 
total transport can be described by the superposition of the diffusive 
transport and a drift component . In his study of flow through long 
tubes, Weber (16) applied the idea of the superposition of a diffusive 
component, a slip contribution, and the viscous component. Weber 
demonstrated that his solution has the 'correct limiting values for 
continuum and free molecule flows and adequately describes Knudsen T s 
data for long tubes (23). Milligan experimentally verified, with 
excellent agreement, the analysis technique of Weber for rarefied flow 
in. long tubes (17) , Lund and Berman .(24) developed empirical relations 
for the flow and self-diffusion of gases in both 'long and short 
capillaries by the superposition , of the diffusive and drift components. 
They developed an .algebraic expression which permits the direct 
computation of the Weber diffusion coefficient at any pressure and thus 
avoided the numerical integration inherent in the Pollard and Present 
treatment of self-diffusion. Lund and Berman demonstrated the adequacy 
of their model for describing self-diffusion and flow in capillaries 
in addition to the flow between flat plates. Milligan, Cowling, and 
Wilkerson (14) extended the superposition analysis technique of Weber 
to long annuli with continuing success. 

In this section, the self diffusion flow in both the grooves and 
annular space of the viscoseal will be discussed and evaluated using 
the Pollard and Present treatment as applied by Weber. 
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Annular Space Self-Diffusion 

The self-diffusion flow in the clearance space of the rarefied- 
gas viscoseal was obtained by considering this flow to be that' of a 
concentric annulus. The molecular transfer was determined by evaluating 
the net number of molecules crossing a plane normal. to the annulus. 

The evaluation was done by considering separately the molecules t which 
come from the outer and the inner, walls with the grooved inner wall 
taken as being smooth and of diameter D as shown in Figure 1. The 
details of this derivation are presented in Reference (3) and the 
procedure used was essentially that given by Weber (16). The following 
equations were numerically integrated to obtain the diffusion 
contribution in the annular space: 


N 


net outer wall 


, ^ R„ 0 ir/2 

= / /°/ <*>• 


R 1 0 0 


sin 0 cos 2 0(1 - e” R/, V dR' d0 d© 


(58) 


N 


net inner .wall 


, _ R„ it tt/2 

- - 2?A |r< r 2 - r i '> ! !.! < E '>- 


E i t 0 


sin 0' cos 2 9'(1 - e dR' del de'. (59) 


The total self-diffusion flow in the annular space is the sum of 

Equations '(58 and 59) and may be arranged into specific molecular 

flow rate form using 

dn _ d__ = 1__ dP £ 1 __ _£. 

dy dy KT dy KT L 


and 

V m 2 r— 1 1/2 
L TTm J 
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to obtain 


N 

AP 


= - 40 


2 ^ 1/2 


( r 2 “ r l 5 ' 


TTinKT 


N k [SUM (O.W.) + -SUM (I.W.-)]., 


(60) 


where SUM (O.W.) and SUM (I.W.) designates the .numerical integration 
of the integrals of. Equations (58 and 59), respectively. 

Groove Self-Pi f fusion 

The self-diffusion flow in the grooves of the shaft was obtained 
in a similar .manner to that of the diffusive flow in the annular space. 
The flow in the groove was determined by considering this transport to 
be that of a long .groove 'of rectangular cross section. This neglects, 
any curvature effects of the helical groove and becomes increasingly 
in error as the groove dimensions become of the order of magnitude of 
the seal diameter. The development work of the numerical scheme is 
presented in complete detail in Reference (25) . 

The' groove diffusion flow involves the 'numerical solution of the 
following type of equation: 


N 


TO , h b’ h 

Kg dn t f f f r (b -rQb 1 

" “ 2ir > * > * s 


0 0 0 0 


(1 - e R/N Kg) 


+ S-. (i - e RR/N Kg)} £ 2 da f (drjdzdQ , 
rr 


where 


(61) 


N Kg b f 

is the Knudsen number of the .groove 'flow with the groove width* b’, 
being taken as the characteristic, dimension. Again 

V = 2[^SS.] 1,/2 
Tim 

and 
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dn = l./'L-'i = ap 

d£ d^Kr KT d£ 'VKT i * 

where SL is now the groove length. The relation between the length of 
one groove and the axial seal length is 


£ _1_TO L L 

Ng cos a (w + b) sin a * 

where N c is the number of grooves. For a multiple grooved shaft the 
o 

total groove diffusion flow is Ng times the flow of a, single groove. 
Thus the total specific molecular flow due to the 'diffusion is 


N 

2P = " 



sin ct 

( — £_ 

7TL TTmKT 



[SUM (Walls)] 9 


(62) 


where SUM (Walls), designates the numerical integration of the integral 
terms of Equation (61) to include the total walls of the grooves. 

2.25 Composite Solution . The total specific molecular flow for 
the composite solution is obtained by adding the continuum solution 
plus the slip flow contribution after correction for the molecules 
which do not experience .slip plus the self-diffusion flows in the 
annular space and the grooves. Thus for the mean-value solution 


composite ‘ ^ < 57 > + E< >' 


(60) + Eq. (62) 


(63) 


It should- be noted in simulating the flows that the contributions of each 
equation must be evaluated at the .same physical gaseous state. In ■ 
the development of the groove self-diffusion the Knudsen number, , 
was based on the groove 'width , b 1 , rather than the .seal radial clearance, 
c. Thus the corresponding value of groove Knudsen number for a given 
clearance Knudsen number is 

H Kg - Vfrf • 


(64) 
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This same technique can be used to formulate the composite 
solution for the variable mean free path analysis and other models 
such as the annulus-groove model,, 

2.3 Annulus -Groove Model . In this section a simpler model than, 
the Modified-Bodgson Model will .be developed which will lend itself to 
a less arduous solution in the transition regime and will also 
correctly predict the continuum performance. The' model also predicts 
the rotor induced flow in. a more rigorous manner. 

The basic model will be essentially the same as the one used by 
Hodgson and which is basically the same as the simplified screw extruder 
theory presented by Carley, t et al. (26). The assumptions inherent in 
these models are: 

(1) the total flow in the seal can be treated as, the super? 
position of the leakage flow in the grooves,, the leakage 
flow over the lands and the rotor induced flow in, the 
grooves , 

(2) the pressure varies continuously along the groove and 

is ‘constant over the cross section of a .particular groove s 

(3) the groove depth is small compared to the diameter of 
the seal, thus allowing curvature effects to be neglected 
in the groove , flow .development, 

and 

(4) the flow in a seal with a grooved housing is identical to 
the flow in a seal with a grooved shaft (13). 

Assumption 1 obviously neglects ‘the convective coupling of the flow 
components. The solution without this assumption is extremely complir 
cated for even purely continuum flow (27) . Its exclusion would 
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certainly lead to an even more complex analysis when non-continuum 
boundary conditions are applied. Since the objective is the development 
of a simplified theory, the inclusion of assumption 1 is a necessity. 

One important assumption that is usually made which is not made here is 
that the groove sidewall effects are negligible. This assumption is 
one of the prime distinctions between this analysis and the Slip- 
Modified. Boon and Tal Model in , which this assumption is made. 

The assumptions governing the development of the component 
flows ares 

(1) the flow is steady, constant viscosity, fully developed, 
isothermal and Newtonian with negligible body forces, 

(2) the Navier-Stokes equations with non-continuum boundary 
conditions are applicable, 

and 

(3) the non-continuum boundary conditions can be expressed as (19) 

v| -- . 

‘wall 3n'wall 

In the development of the flow components, flow models which have been 
experimentally verified will be used and' reference made to their 
verification. 

2.31 Flow in . the Groove 

Although most previous investigators .have initially treated the 
groove flow in two parts, this analysis will initially treat the groove 
flow as a single flow from which the two previously mentioned groove 
flows are eventually obtained. 

Based on the assumptions stated above, the Navier-Scokes equations 
reduce to the. single z-momentum .equation 



-32- 


3 2 v 3fv = JL dP 

3x 2 3^ 2 y dz * 

The applicable slip boundary conditions .shown in Figure 5 are: 

% co.7> " 0 


v(b’/2,y) = - GA — (b*/2,y) 

3x 

v(x, 0) = GA ~ (x, 0) 


_ 

v(x, h) - - U cos a - GA-^Cx, h) 


where advantage -has been taken of the symmetry about the y-axisi The 
slip coefficient* G, is' usually taken as unity as will be done here. 

It is convenient at this point ,to non- dimens i on alize the velocity 
and the coordinates and introduce an index of rarefication. The non- 


dimensionalized variables are taken as 


x = x/ (b * /2) 


y = y/h and 


U cos ct* 


A common index of rarefication, the Knudsen number is taken as 
N Kh = ^ h * 

Introducing these new variable into Equations (65) and (66) one 


obtains 


3 2 u .a. 2 £jx _ (b T ) 2 dP = Rj_ 

3x 2 3y 2 ’4]iU cos a dz U 


(b*) 2 dP 

4ycos a dz 


and the , boundary conditions are 


S <°» = 0 

n > 2N Kh 8u n \ 
u(l, y) — (l, 7) 
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u(x, 0) = N Rh (x, 0) (68c) 

and 

u(x, 1) = - 1 - (x, 1). (68d) 

Many techniques have been used to solve non-homogeneous problems 
of the type ■ presented above. The- least arduous is a modified form of 
the method of variation of parameters (28) . This method has been 
employed successfully by Ebert and Sparrow (29) to solve 1 Equation (65) 
with' four homogeneous boundary conditions. In -the case considered 
here Equation .(68d) presents a non-homogeneous boundary condition,, but 
the method used by Ebert and Sparrow still leads to a solution. , 

The method of variation ,of parameters is a generalization of the 
method of separation of variables, and as, such, experience gained from 
employing the latter is helpful when applying the former method. In 
the method of separation of variables, the sign of the .separation 
constant is taken so that the trigonometric solution is obtained for 
the homogeneous direction. In 'the method of variation of parameters 
the same reasoning is used in the assumption that the solution has the 
form 

CO 

u(x, y) = Z 0_(y)(cos a x + C sin, a x) (69) 

1 ** n / n , 

where the x-direction is the .homogeneous direction determined by ,an 

inspection , of Equations (68). The function 0 (y) is a yet unknown 

n 

function of y. One now proceeds to determine 0 „ a B and C_. 

n® n®, 7 

Substituting Equation (69) into the boundary condition of Equation 


(68a) one obtains 
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03 

2 0 I1 (-y)[ c 7 a n cos (°) - a n sin (°)] = Q 

n = 1 

which implies .that.C^ = 0 and that 

CO 

u(x, y) =' 2 0 (y) cos a x.. (70) 

^ n / xi 

n = 1 

Substituting Equation (70) into Equation (68b) gives 

® 2N a 

2 0 (y) [ cos a — sin a ] = 0 

n a n 

n = 1 

from which it follows that for a non-trivial .solution . 

2N a 

cos a — sin a = 0 

n a n 

or 


a tan a 
n n 



(71) 


The eigenvalues, of the eigenfunctions, cos Ct^x, are the roots of 
the transcendental relationship of Equation (71). 

The problem now reduces to determining 0 n (y). Substituting 
Equation (70) into Equation (67) yields , 


CO CO 

2 0 ^ 0 (y) cos a x + (— ) 2 2 0"(y) cos a x 

-nn n . / _n n‘ 

n = 1 n = 1 


= 2 cos a x 

TJ _ n ' n 

n = 1 


(72) 


where satisfies the .Fourier series 
n - 


1 = 2 cos a K. 

n n 

n = 1 


(73) 


The Fourier cosine series coefficient, must, be 

2 sin a sin a 

q _ 2 _ 2_ r n , 

n„ a + sin a cos a , a L ' 7 2K, sin^ a 
n n- n nl+h n 


(74) 


a 
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where the last step comes from Equation (71), 


In order for 0^(y) to satisfy Equation (72), it must be that 
2a 4R’ft 

0;(y) - (-^) 2 0 n (y) - —T * 

Ua 


( 75 ) 


The determination ,of 0^ thus reduces to the solution of a non- 
homogeneous second-order ordinary differential equation,, The boundary 
conditions on 0^(y) are obtained from Equations (68c) and (68d)„ * 


From Equation (68c) , one boundary condition is 

.0 (0) = N0'(O). 
n Kh n 

Substituting Equation .(73) for the unity term in Equation (68d) 
leads to the second boundary condition. 


(76a) 


* - S n - 


(76b) 


The solution of Equation (75) subject to the boundary conditions 
in Equations (76) is 


-9 y -9 y R'fl 

0 (y) «-(C, + C. )e n + (C„ + C_ )e n ~ 

v n vj/ v lnp lnu 2np 2nu ^2 


where 


2a 

9 n= a 


n 


-9„ 


r R ' a n ■ - - U + Vn) 

lnp „_.2 „ v2 -0 n ~ n2^© 


n (1 - W e - (1 + \h e n> e n 


n 


(77) 


a (i + 0 ) 

n Kh n 


lnu 


J 2np 


«• - V,) 2 « A - < x + w 2 * 8 ” 

r d 1 ”. 

(1 - N , 9 ) 2 e“ Sn -- (1 + 'N 0 ) 2 e® n Ua 2 
' Kh n ; ' Kh n n 
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o (i - e ) 

Kh n 


"2nu 


(i - N^ejV 6 * - (i + N Kh e n ) 


2 e. 


Kh n 


n 


By substituting Equation (77) back into Equation (70) the 
dimensionless velocity distribution becomes, after considerable 
simplification , 


u(x, y) = 


(b') 2 


dP 


CO 

E [- 


sin a 


n , 


2rn cos a ds n ; !*! + 2%.,^ „ 


]. 


COS 0t X 
[ 1 ]]• 


2a v a 

. n y , , n . , 2a y 

cosh tanh — smh n -7 

a a — ■ — 


or 


n 


1 + 


2N K& 


tanh 


■^T 


- i] 


CO 

- E 


sin a 


2 [- 


n 


2N 2 

n - 1 1 + Kh sin a 


cos a x 

JE— — -] 


a 


n 


2a y 2N,„_a 2a y 

n J Kh n , 

smh 1 ■ - + — cosh ■ — 

a a a 


n 


2a 


r 2 N_ a 2a 4H_.a 

r , . ( Kh n v 2, , , n , Kh n • n 

[1 + v — ) ] sinh + cosh. 

a a a a 




(78) 


The volume flow rate in the groove is obtained by integrating the 
velocity distribution over the area of the groove. The volume flow 
rate thus obtained is 


v , = (b * ) 3 h dP ; a 
6 2y dz n = 1 a 5 


sin 2 a 


E- 


n 


2N T1 


-][■ 


tanh a /a 
n 


• , , Kh . 2 .. . . 2N Kh K n , . n 

n 1 + r- sm a 1 + tanh — 

a " n, a a 


a 

E] 

a a J 


• 2 „ 
sm a 


= Ub T h -cosa 


n 


n = 1 a~ 


][■ 


tanh a /a 

n 


n 


1 + 


2N n J L 2N_ a 

Kh , 2 -t , Kh’ , , n 

sin a 1 + tanh — — 

ana a 


- 3 . (79) 


The first summation term represents the pressure induced flow in the 
groove whereas the second term represents the rotor induced flow. 
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The linearity of Equation (65) has .been demonstrated in the, 
analysis above since the 'solution of Equation (79) can be shown to be 
the sum of two solutions." - The first term is in agreement with the 
expression obtained by Ebert and Sparrow (29) for slip flow in a 
rectangular -duct with stationary walls. Milligan and Patterson .(5) 
have experimentally verified the solution of Ebert and Sparrow. The 
last term. is the solution of V v “ 0 with the boundary conditions , given 
in Equation (66). It is thus demonstrated that it .is permissible to 
obtain the two solutions mentioned above -independently and then to 
add them together to arrive at the same result as Equation , (79) „ 

At this .point it is interesting ,to determine V 1 when the Knudsen 

O 

number approaches .zero, the .continuum flow regime. The solution for 

the eigenvalues becomes r- 

cos a = 0 
n 

or 

« n = 2|‘ n = 1, 3j 5, 7 . . . 


With these eigenvalues and with ;N Tr , 0. V r becomes 

Kh g- 


V 


o CO 

, 16 (b f ) h dP v a r . , mr mr, 

* = 3 -T- S —j[tanh "Jr “ 2^3 

g 7 r\l dz n = 1. n Za 


8Ub'h cos a 


IT 


v a . , mr 

_ _ , n 3 nh 2a 
n = l n 


n = 1, 3, 5, 7 


0*0 


(80) 


2.32. Correction ,to Rotor Induced Flow . In the analysis above 
it was assumed that the rotor velocity, U cos a, acted at the 
top of the groove (y = h) . No attempt was made to account for 
the fact that the .rotor is not -located at y = h but rather at 
y = h + c. Since an exact analysis of the region .above the 
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groove, h_<y<_h + c » has been shown to be very complex 'even for 
continuum flow '(27) , some sort of approximation must ’be .made to obtain 
a simple solution. The following assumptions are made: (1) the flow 
induced by the rotor in 'the region b <_ y <_ h + c does not contribute 
to the seal discharge and, (2) an effective inducing velocity less than 
U cos a acts at y = h as a result 'of the clearance region h <_ y <,h + c. 
This effective inducing .velocity U e , will be determined by calculating 
the average velocity in the plane y = h of a, hypothetical groove of. ' 
depth h + c.and width b with the upper boundary moving at U cos a and 
with slip boundary conditions on all surfaces. The hypothetical' 
groove and associated boundary 'conditions are shown in Figure 6. The 
brackets in Figure 6 indicate the portion of the boundary over which 
each boundary condition is, applied. 

From Equation (78) the velocity distribution 'induced by the, rotor is 


sin oi 


E — [- 

i a 

n ,= 1 n . 


cos a x 
n 


2N Kg'n, 


2N a 2a 4N o; 

) 2 ]sinh —a +T J SL* 
a a a 


If Equation (81) is .applied to a 'groove of depth h + c, the resulting 
distribution is 



[1 + (■ 
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where 


a' = 


, f 1 


h + c 1 + c/h 


N. 


= 


Kh 


Kh f h + c 1 + c/h 


y » = JL . 

y h + c * 


The eigenvalues of Equation (71) are still the same because 


2N. 


Kh 


21 W * 


The average velocity at the plane y = h or y* - 


U e - / v^(x,h/(h + c))dx. 


Carrying through the integration, U e becomes 


sin 2 a 


U = 2U cos a 2 -~[- 

e n - 1 a 2 


n 


2N 


. - . Kh . 2 „ 

n 1 + sm * a 

a n 

- 2a 2N t _, a 2a 

. , n . Kh n , n 

sinh — — + — cosh — — 

a a a 

2N a 2a 4La 2a 

[1 + — ) 3 sinh r + - cosh ~~r 

a a a a 


or 


U e = [U cos a] £ c 


where the 2 c is equal to twice the summation in Equation <(82) 
corrected groove flow is now obtained by replacing U cos ct in 
(79) with [U cos a]E c « The total groove flow now becomes 


V' 

g 


( b T ) 3 h dP 
2U dz 


Z - U b'h cos a[Z 2 ] 
p c u J 


(82) 

(83) 
The • 

Equation 

(84) 


where 
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2 = 

P 


CO 

2 

T1 — 1 


sin 2 a 


[" 


n 


H- 


a 

tanh — 
a 


a 5 \ . 2N Kh . 2 .. J \ . 2N KhV. , a n 

n i + — — sm a 1 + tanh — 

ana a 


a 

ti 


] (85a) 


CO 

2 


sin 2 a 


a 

tanh n 


S “ = nil a 3 ' + 


2W. 


-][_ 


Kh . 2 „ 2N V , a a 

- — sin , a - , Kh n . . n 

a n 1 + — — — r- tanh - — 

a a 


(85b) 


and 


CO 

2 = 2 
C n = 1 a 2 


. 2 C4 
sm n 


2N, 


J. 


n 1 + 


Kh . 2 „ 
— sin a 
a n 


2a 2N T „ a 2 a 

. , n , Kh n , n 

sinh — — + — — cosh 

a a a 


2N_ , a „ 2 a 4N ,a 2a 

[1 ( _ SL n )2] h _u + _SLn _n 

L v a J a. a a' 


)• 


(85c) 


In terms of the viscoseal geometry 


dP 

dz 


dP 

d"£ 

g 


N g (b t + w) 
1TD 


dP . „ dP 

dL = sln “ dL 


( 86 ) 


and since the pressure gradient in the viscoseal is assumed to be 
constant , 


dP 

dz 


sin 



(87) 


The total groove flow is T obtained by multiplying Equation (84) by N k . 

s 

Noting that N = (TD) sin a/(b’ + w) the total groove flow becomes 

tt _ „ tti _ Trl^b 1 ) 3 bAP sin 2 a v UfrDb'h cosa sina r r 

V g " N s V g 2 vTl (b * + w) L p (b r + — u o’ <88) 

2.3 Land Leakage Flow 

The land leakage flow is taken* as flow through a long annulus. 
Milligan et ,al. (4) have developed an expression for slip flow in an 
annulus and have obtained excellent experimental confirmation of the ‘ 
expression., The expression derived by Milligan is 
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7 - - , 1,1)4 . i%- A, 

L 128}i d& L L 

where 

*L ° Itl - K'h + + 

[2K t (K ? 2 - 1). - 2K 1 (In K’)(k' 2 - 
+ 2(1 - K' 2 )(l k') 2 n k ], 


2N k (1 - K' 2 ) 

K' lu K' - N (1 - K’ 2 ) 
Jx 


K T + 1) - 


(1 - K >2 ) 2 
2 In K T 


( 89 ) 


(90)- 


and is .the effective land pressure gradient. For a single-threaded 

Li 

seal, Hodgson (9) showed based on , assumption (2) at the beinning of 
this chapter that 



w’ + b* 
w' 


cos a 


dP 
dL * 


(91) 


As was shown in Reference (13) , Equation ,(91) remains valid -for a 
multi- threaded seal. 


Combining Equations (89) and (91) the land leakage flow becomes 




TrD^(w' 


+ b*) cos a 


128y w 1 




dP _ 7TD^(w 1 + b 1 ) cos a 

dL “ 12 8p w'L 


*L 


AP 


(92) 


2.3 Total Seal Flow and Sealing Coefficient 

The total seal flow ds obtained by combining Equations (88) and 
(92). The- result is 


AP 

V = [A E 
T liL p p 


+ A 


: V - 


6U 

2 


[ASS] 
u .c u 


(93) 


, 7rD(b’) h sin z a 

p = 2 (w ' + b 1 ) 


where 
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TTD 4 (b t + W*) COS a 

128 w’ 


and 


A u 


'rrpb'he^ co9 a sin a 
6(w’~'+ b') 


The specific molecular flow rate is 

1L = Zl C JE) 

Ap Ap v, kT'' 


2iP “ ULkT ^ A p S p + A c A L‘* “ 2 t ^ - lU A u S c 2 J 


(94) 


or 

r + 1 , 

[-£ 
c^kT r p 

where Equation (12) is >used to eliminate Ap„ 

A common- dimensionless viscoseal performance index is called the 
sealing coefficient and is defined by 


A = 6ym 

c 2 Ap (95) 

where AP is the pressure difference at zero flow., By equating to 
zero in Equation (93) , the sealing coefficient is obtained as 


A 


A 2 + A A 

p p c L 

"'AH 


(96) 


u c u 

A digital computer was used to compute the specific .molecular 
flow rate from Equation (94) and ,the sealing coefficient from Equation 
(96). In order to evaluate the summations. S s S , and I . it is 

p» c » u 

necessary , to determine the eigenvalues which satisfy Equation (71). 
Newton’s method of approximating roots ,was employed to solve Equation 
(71) for these eigenvalues at eqch Knudsen number. 
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3. EXPERIMENTAL INVESTIGATION 

Since no reliable performance data were available for visco.seals 
having a rarefied gas as the sealant it was necessary to conduct an 
experimental investigation. The object was to obtain performance data 
which could be used to evaluate the theoretical models . These 
investigations were conducted on multiple groove two-inch diameter 
viscoseals operating over a wide range of shaft speeds and gas densities. 
Data were obtained with no leakage through the seal which permitted the 
evaluation of sealing coefficients. In addition data were obtained for 
a wide range of net leakage conditions . All rarefied data were obtained 
using argon as the sealant. Continuum sealing coefficient data were 
obtained using air as the sealant. 

3.1 Viscoseal Test Section 

The experimental apparatus was designed to investigate viScoseal 
performance in the gas flow regime between continuum and free molecule flow. 

The viscoseal test section, Figure 7, consists of an outer housing 
with its. associated vacuum pumping system surrounding a rotating grooved 
shaft. The shaft is a hollow eight-inch cantilever extension of a high 
speed spindle shaft which is belt-driven through an intermediate spindle 
by a direct current motor. The drive system is capable of seal shaft 
speeds from zero to 35,000 rpm. The speed control for the motor is 
self-regulating and maintains a selected speed within +0.1 percent. 

Since the entire test section operates under vacuum, a rubbing contact 
graphite ring seal is provided where the rotating shaft penetrates the 
housing. The ring seal is a series B-103032, Type E, manufactured by 
the Cleveland Graphite Bronze Division of the Clevite Corporation, 

Cooling water and air are supplied at the ring seal end of the spindle 
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to minimize thermal growth, due to the rubbing seal friction ; under 
dynamic conditions. A Gonflat vacuum flange is provided between the 
housing and its support for purposes of sealing and two provide minor 
adjustment capabilities for housing-shaft alignment. Other vacuum 
seals are obtained by the use of "0" rings in addition to vacuum 
sealants for threaded connections. The maintenance of 'a high quality 
vacu um system with essentially no atmospheric leakage, except at -the 
shaft seal, was assured by the frequent "use of a helium leak detector 
throughout ,, the experimental program. 

Experimental data were obtained for two different viscoseal 
geometries consisting of a grooved shaft inside a smooth housing. 
Pertinent -specifications of the rarefied viscoseals are contained in 
Table I, 

Table I. Dimensional Specifications of Viscoseals (All Dimensions s/re in 
Inches) 


Parameter 

Seal No. 

Housing diameter (+,0.0003) 

2.0088 

Shaft diameter (+ 0.0001) 

2.0005 

"Gold" radial clearance 

\ 

(+ 0.0002) 

0.00418 

Axial length (+ 0.005) 

4.530 


1 S eal No. 2 

2.0088 
2.0005 

0.00418 
4.530 

0.1691 + 0.012 
0.2016 + 0.012 
0.0144 + 0.0005 
19.474° 


Groove axial width 
Land, axial width 
Groove depth 
Groove helix angle 
Aspect ratio, (b cos a) /h 


0.03111 + 0.0003 
0.03235 + 0.0003 
0.03065 + 0.0003 
9.30° . 

1.002 


11.071 
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A schematic diagram of the overall experimental test apparatus 
is shown in Figure 8. 

3.2 Vacuum Pumping System 

It was necessary to develop .a vacuum pumping system capable of 

providing the desired low pressures and gas flow rates. This task was 

accomplished by connecting two independent vacuum pumping systems in 

parallel to a common , reservoir. Figure 8. One of these pumping systems 

consisted of a single stage rotary oil-sealed mechanical pump mated with 

a three-stage high vacuum oil diffusion pump. This part of the pumping 

system was connected to the reservoir with a large diameter open-or-close 

vacuum valve. The second pumping system was composed of another single 

stage rotary oil-sealed mechanical pump mated with a positive displacement 

roots blower type vacuum pump. This portion of the pumping system was 

connected to the reservoir ,througha throttlable high vacuum valve. 

Through ^manipulation of the too connecting valves it was possible to 

regulate the pumping speed over a suitable range of downstream test section 

pressures. The complete vacuum pumping system was - capable of attaining 
-5 

pressures to 10 millimeters of mercury with a blanked-off system. 

3.3 Ins trumentation 

The radial alignment of the shaft within the housing was determined 
by the use of five proximity detectors manufactured by the Bentley- 
Nevada Corporation. These probes were located near either end of 
the viscoseal section. When viewed from the shaft drive end of the 
test section, the probes of the first set were located at 12, 3, 
and 6 o'clock positions, and the probes of the second set were 
located at 12 and 3 o'clock positions. The arrangement of the 
probes permitted the shaft to be aligned within the housing to a 
value of eccentricity ratio (shaft centerline deviation/mean .radial 



-46- 


clearance) on the order of 0.05. The. probes were also used to determine 
the thermal growth of the shaft caused, by the graphite ring seal 
friction when the shaft "was rotating. In addition,, the probes 
permitted the measurement of • all vibrational movements of the shaft. 

The output gains of the detector amplifiers were individually adjusted 
and calibrated to ensure a linear output voltage of the probes as a 
function of the clearance gap. 

Pressure measurements upstream and downstream of the viscoseal 
were obtained using both McLeod gauges and absolute aneroid type 
gauges. The McLeod gauge is normally considered to be a primary -standard 
(30), but 'after previous attempts (31) to use thermocouple gauges, 
ionization gauges and cold cathode gauges, it was apparent that McLeod 
gauges were the .only instruments capable of giving the desired accuracy 
and reproducibility. The McLeod gauges utilized were the, GM-100A 
gauge manufactured by the .Consolidated Vacuum Corporation and are 
described in detail in -Reference (31). 

Temperature measurements of the argon entering the high pressure, 
end of the viscoseal test section .were made early in the experimental 
program using a thermocouple, A comparison of the measured gas 
temperature with the ambient (room) temperature revealed that -the error 
involved in using room temperature instead of the actual gas temperature 
was less than one percent. Therefore the gas temperature was taken to 
be 297°K in all data calculations. 

Speed measurements of the viscoseal shaft were made using a 
magnetic -pickup located near the attachment nut for the drive belt 
pulley on the high speed spindle. The pulses from the pickup were 
registered using an electronic counter to obtain shaft rpm. The speed 
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measurements were , independently verified using a calibrated hand held 
tachometer in. the low rpm range and a calibrated strobotac over the 
entire speed range of 0 to 35,000 rpm. 

The molecular flow rate .through the viscoseal test section was 
determined using a variation of a constant pressure method developed 
by J. R, Downing (30) for the measurement of pumping speeds of vacuum 
pumps. The technique consists of applying the perfect gas equation 
of state to known gas volumes at two different .time periods. In brief, 
the flow measurement is .obtained by adjusting the indexed valve of 
Figure 8 to a desired setting and waiting for steady state to be 
achieved in the test section .as indication .by pressure and proximity 
probe measurements. The shut-off valve leading from the constant 
pressure gas reservoir is then closed and the .rate of rise of manometer 
fluid in the right hand side of the manometer is observed. With 
knowledge of the manometer cross ' sectional area, the rate of fluid 
rise, the value of the calibrated volume, the initial and final pressures 
in the calibrated volume, and the gas temperature, the molecular flow 
rate is calculated in the following manner. 

Consider the flow measurement diagram of Figure 9. Let 

V q = calibrated volume including tank and manometer down to 
the zero deflection line, "0" 

P q = atmospheric pressure 

= initial manometer deflection , 

= final manometer deflection 

A„ = cross section area of manometer tube 
m 

W = specific weight ofmanometer fluid 
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From fluid statics:- 


P, = P + 2WH- 

1 o 1 

P~ = P„ - 2WH 

2 o z 

The corresponding gas volumes are; 


“ V 

+ H 'A 

o 

1 m 

= V 

- H„A 

0 

2 m 


Applying the perfect gas equation of -state 


5 


PV = NKT, ( 97 ) 

on a molecular basis to determine the number of molecules within the 
system ,at time ’( 1 ) and time ( 2 ) , then. 


4H " N 1 - » 2 


P V 
11 

KT, 


V2 

KT, 


(98) 


1 iVi 2 

Assuming isothermal conditions and substituting for the pressures and 
volumes. Equation (98) simplifies to. 


AN = 


H l + H 2 

KT 


{ A m [P o + 2W(H 1( - H^) ] +2WV } 


(99) 


The flow measurement system was operated experimentally such that* 

H 1 = H 2" 

Letting 

H = H ! + H 2 

and At.be the time period for the total deflection ;H, then the molecular 
flow rate is 


w At 


H 


AtKT 


{A P + 2WV }. 
mo o 


( 100 ) 


The basic premise for the .measurement system rests in. the fact 
that while the pressure in, the calibrated volume, does change slightly, 
less than + 1 %, it -still remains very close to, atmospheric pressure. 
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Thus, the change in flow rate through the valve is negligibly affected 
by the slight decrease in the pressure within the calibrated volume 
during the course of a measurement. 

3.4 Experimental Procedure 

Experimental data of two different types were obtained for the 
operation of the viscoseals in the rarefied regime. These may be 
described as net flow leakage data and sealing coefficient data. The 
operational procedure utilized in gathering the data differed only in 
respect to shutting off the argon gas supply when obtaining sealing 
coefficient data. 

A.11 data were acquired by the coordinated adjustment of the gas 
supply rate, the vacuum pumping speed, and the seal shaft drive speed 
followed by a time period sufficient to obtain a steady-state condition. 
Following achievement of steady state, final readings of the gas flow 
rate, pressures upstream and downstream of the viscoseal, and the 
proximity probes were made. The proximity probes were used to record 
the seal clearance which decreases slightly with increasing shaft 
speed due to the friction of the graphite ring vacuum seal. 

3.5 Data Treatment 

The viscoseal performance data for both the net flow leakage tests 
and the sealing coefficient tests were correlated versus an index of the 
flow rarefication, Knudsen number. The mean radial clearance, c, was 
selected as the viscoseal characteristic dimension. Thus using the 
Chapman relationship for the mean free path, the experimental values 
for Knudsen number were determined using the following equation. 

R 

x Vtt* ir T>1/2 


C 


2P c 


( 101 ) 
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The mean free path was based on. the average seal section pressure 


P = 


P T +P b 

2 


( 102 ) 


The molecular flow rates for the net .flow leakage data were 
calculated using Equation (100). These data were revised to specific 
molecular flow rates .by dividing by the pressure drop, AP, across the 
viscoseal, or 


' H[A P + 2WV ] 
N mo o 

AP = (P T - P b )AtKT 


(103) 


The sealing coefficient data were calculated from the zero net 
leakage tests using 


A = Mli = 6JJUL . - (104) 

APc (P T - P b ) c 

3.6 Experimental Uncertainty . A detailed analysis was performed 
to estimate the propagated uncertainty in the final performance parameters 
based on the individual uncertainties of all measured variables. The 
analysis was performed using standard statistical techniques after 
postulating a normal distribution in the uncertainty of each variable.' 
This method of analysis assumes all errors are random in nature. 
Considerable effort was expended during the experimental investigation , 
to minimize -any error of a systematic nature. 

The propagated uncertainties, in general, are not constant and vary 

with the .degree of rarefication., The uncertainty in the specific 

* 

molecular flow rate, N/AP, of Equation , (103) is + 5.1“percent at an 

Inverse Knudsen Number of 27 where the calibrated volume, V , 

o’ 

contribution is 84 percent of the total. The uncertainty in the. 
specific molecular flow rate increases to + 8.9 percent at an Inverse 
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Knudsen Number of 0.50 where the uncertainty in the differential 
pressure, AP, contributes , 72 percent of the total. The ability to 
determine the Inverse Knudsen Number of Equation (101) is estimated 

at + 4.8 percent for 1/N equal to 27 and increases to 4- 5.7 percent for 
1/N equal to 0,50. The propagated uncertainty estimate in the sealing 
coefficient, A, of Equation (104) is insensitive to the degree of 
rarefication and remains at essentially +9.6 percent. The uncertainty 
in the radial clerance,, c, contributes -approximately 99 .percent , of 
the propagated uncertainty estimate for A. It should be noted that 
the clearance uncertainty is controlled by the ability to establish 
the diameters of the seal shaft and housing and not by any limitations 
of the proximity probe system. 

All uncertainties are stated to a confidence level of 95 percent. 

3 . 7 Investigations Conducted 

Sealing coefficient performance data and net leakage performance 
data were obtained for each of the two geometries . Sealing coefficient 
data were obtained over a wide range of sealant gas rarefication for 
shaft speeds of zero, 5000 „ 10,000, and 30,000 rpm. Net leakage data 
were also obtained over comparable ranges of rarefication and shaft 
speeds. These data are presented and discussed in Section 4 and the 
reduced data are tabulated in Appendix 'B. 

4. RESULTS 

Two types of results will be presented for each of the three 
analytical models.' Net leakage and sealing coefficients will be given 
as a function of the degree of rarefication. The inverse Knudsen 
number is used as the parameter indicating the degree of rarefication. . 
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4,1 Net Leakage 

-The net leakage -f!low occurs from the high pressure end to the 
low pressure end of the viscoseal when the viscous pumping 
action is insufficient -to overcome the pressure- induced flow in the 
grooves and through the annular clearance space. The net leakage .data 
are shown k In -Figures 10 through- 15. These data are presented in the 
form of specific leakage rate .(units of molecules per sec-yHg) versus 
Inverse Knudsen Number. The experimental data 'for viscoseal Seal No. 1 
are shown 'in Figures 10, 11 , and 12. The theoretical specific leakage 
values as. predicted by the Modified Hodgson Model are show in 
Figure 10. These theoretical values were obtained by using Equation (6) 
for the zero rpm case and Equation (8) for the 5000; 10,000; and 30,000 
rpm cases. The Slip Modified Boon and Tal Model was used to determine 
the theoretical, values show on Figure 11. Two solutions for this model 
are presented. The mean-value values were calcualted using Equation 
(55). The variable-mean- free-path values were obtained by solving 
Equation (46). In Figure 12 theoretical values for the Annulus-Groove 
Model are presented. Equation (94) was. used to -determine the theoretical 
values show in Figure 12. 

The experimental data for Viscoseal No, 2 are presented in, Figures 
13, 14, and 15. These .figures contain data for one , additional speed, 

2500 rpm. The theoretical specific leakage values show. on these 
figures were determined. in the same manner as described above for 
Viscoseal No. 1. 

In Section 2.25 a composite type solution was discussed for the 
Slip Modified Boon and Tal Model, This type of composite solution ,is 
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obtained by the superposition of the corrected continuum-slip flow 
contribution with the self-diffusion f lows , in the annular , space and 
the grooves. The corrected continuum-slip flow is determined by 
taking into account the fact that all the molecules cannot , experience 
the slip boundary condition. This procedure is explained in detail 
in Section 2,23. A detailed procedure for calculating the self-rdif fusion 
flows is given in Section 2.24. Tables II through VII present composite 
solution results and the uncorrected continuum-slip results. In each 
table the first column on the left is Inverse Knuds en Number, Moving 
to the right in the tables the second column is the uncorrected 
continuum-slip solution., The third column is the , particle-slip 
corrected solution. The fourth column is the self-diffusion through the 
annular clearance space, the fifth column is the self-diffusion through 
the , grooves , and the last column is the composite solution. 

Table II shows specific leakage results from the Slip-Modified 
Boon and Tal Model for Seal No. 1 at - zero rpm. Table III 'is 'the same 
except that these results are for 5000 rpm. Tables IV and V contain 
specific leakage as predicted by the Boon and Tal Model for Seal No. 2 
at zero and 2500 rpm, respectively. Tables VI and VII contain results 
for the Annulus -Groove Model, Table VI is for Seal No. 2 at zero rpm 
and Table VII is for the same seal at 2500 rpm. These are typical 
composite solution results. The particular results shown in Tables II . 
through VII were chosen because they were the cases which gave the? 
largest" deviations between the composite solution and the simpler non? 
corrected continuum- slip solution, , 



Table .II. Slip Modified Boon apd Tal Model 

-14 

Net Leakage Solutions x 10 for Seal No. 1' at ’Zero rpm 


ft/AP 


$/AP Composite Solution - Eq. 

(63) 

1/N K Cont. + slip' 

No. Part. Corr. 

Cont. + slip 
Part. Corr, 

Annulus 

Diffusion 

.Groove 

Diffusion 

Composite Flow 
1+2+3 

1 x 10~ 2 

2.22 

0,198 

1.92 

0.48 

2.60 

3 x 10“ 2 

2.24 

0.514 

1.’69 

0.42 

2,62 

7 x 10" 2 

2.30 

0.951 

1.41 

0.34 

2.70 

-i , 

1x10 * 

2.33 

1.18 

1.30 

0.30 

2.78 

3 x 10 -1 

2.58 

1.99 

0.87 

0.17 

3.03 

-1 

7x10 

3.08 

2.78 

0,57 

0.086 

3.44 

1 x 10° 

3.45 

3.23 

0.45 

0.053 

3.73 

3 x 10° 

5.92 

5.84 

0.19 

0.022 

6.05 

7 x 10° 

10.8 

10.8 

6.090 

0.0096 

10,9 

1 x 10 L ' 

14.5 

14.5 

0.064 

0.0068 

14.6 

3 x 10 1 

39.2 

39.2 

0.022 

0.0023 

39.2 

7 x 10 1 

88.4 

88.4 

0.0094 

0.0010 

88.4 

1 x ,10 2 

125.0 

125.0 

0.0065 

0.0007 

125.0 



Table III. Slip Modified Boon and Tal Model 

-14 

Net 'Leakage Solutions x 10 for Seal No. 1 at 5000 rpm 



N/AP 

N/AP Composite 

Solution - Equation (63) 


1/N No 

K 

i* Part. Corr. 

Part Corr. 
1 

Annulus 
Diff. 2 

Groove Diffusion 
3 

Composite Flow 
1 + 2 + 3 

1 x 10" 2 

-8.36 

^10 , 36 

1.55 

0.48 

-8.33 

-9 

3 x 10 

-8.35 

-10.04 

1.45 

0.42 

-8.17 

7 x 10 -2 

-8.31 

- 9.60 

1.22 

0.34 

-8.04 

-1 . 

1 x 10 

-8.28 

- 9.37 

1.11 

0.29 

-7.97- 

-1 

3 x 10 

-8.07 

- 8.61 

0.76 

0,16 

-7.69 

7 x 10" 1 

-7.61 

- 7. §7 

0.49 

0.08 

-7.30 

1 x 10° 

-7.26 

- 7.44 

0.39 

0.059 

-6.99 

2 x 10° 

-6 . 04 

- 6.14 

0.22 

0.031 

-5.89 

3 x 10° 

-4.81 

- 4.87 

0.16 

0.020 

-4.69 

4 x 10^ 

-3.57 

- 3.62 

0.12 

0.015 

-3.49 

5 x 10^ 

-2.33 

- 2.37 

0.10 

0.012 

-2.26 

6 x 10° 

-1.08 

- 1.11 

0.086 

0.010 

-1.01 

7 x 10° 

+0.16 

+ 0.13 

0.075 

0.009 

+0.21 

1 x 10 1 

3.90 

3.88 

0.054 

0.0063 

3.94 

3 x 10 1 

28.84 

28.83 

0.019 

0.0021 

28.85 

7 x 10 1 

78.73 

78.72 

0.0082 

0.00092 

78.73 

1 x 10 2 

116 . 10 

"116.10 

0.0056 

0.00065 

116.11 



Table IV. Slip Modified Boon and Tal Model 

i 

Net Leakage Solutions x 10 for Seal No . 2 at Zero rpm 



N/AP 


N/AP Composite .Solution - 

Equation (63) 

Cont. + slip, 
1/N^ No. Part. Corr. 

Cont. + Slip 
Part. Corr. 

1 

Annulus 

Diffusion 

2 

Groove 

Diffusion 

3 

Composite . 
Flow 

1 + 2 +.3 

-9 

1 x 10 

2.09 

0.110 

1.92 

0.72 ■ 

2.75 

-9 

3 x 10 

2.11 

0.303 

1.69 

0.49 

2.48 

7 x 10“ 2 

2.14 

0.609 

1.41 

*0.31 

2,33 

1 x 10” 1 

2.17 

0.790 

1.30 

0.24 

2.33 

3 x 10" 1 

2.33 

1.518 

0.87 

0.11 • 

2.49 

-1 

7 x 10 

2.66 

2.21 

0.57 

0.05 

2.83 

1 x 10^ 

2.91 

2.57 

0.45 

0.035 

3.05 

3 x 10° 

4.54 

4.41 

0.19' 

0.012 

4.61 

7 x 10° 

7.78 

7-72 

0.09 

0.005 

7.82 

1 x 10 1 

10.21 

10.17 

0.064 

0.0036 

10.24 

3 x 10 1 

26.41 

26.39 

0.022 

0.0012 

26.41 

7 x 10 1 

58.80 

58.80 

0.0094 

0.0005 

58.81 ' 

1 x 10 2 

83.10 

83.09 

0.0065 

0.00036 

83.10 



Table V. Slip Modified Boon and Tal Model 

-14 

Net -Leakage Solutions x 10 for Seal No. 2 at 2500 rpm 



. 


. 1 



"L 


N/AP 

*> 

N/AP Composite Solution - Equation 

(63) 


Cont.+Slip 

Cont. + Slip 

Annulus ; 

Groove 



1/N„ No. 

. Part; Corr. 

Part. Corr. 

Diffusion 

Diffusion 

Composite Flow 



1 

2 

3 

1.+ 2 + 3 


1 x 10" 2 

-2.26 

-4.36 

1.75 

0,71 

-lo9Q 


3 x 10 -2 

-2.25 

-4.15 

1.52 

0.47 

-2.16 


7 x 10 -2 

-2.22 

-3.82 

1.28 

0.29 

-2.25 


1 x lO -1 

-2.20 

-3.64 

1.16 

0.23 

-2.25 


3 x 10" 1 

-2.07 

-2.69 

0.79 

0.10 

-1.80 


--j 

X 

t - 1 

o 

1 

t - 1 

-1.80 

-2.25 

0.51 

0.047 

-1.69 


1 x 10° 

-1.58 

-1.91 

0.41 

0.033 

-1.47 


2 x 10° 

-0.829 

-1.00 

0.23 

0.016 

-0.75 


3 x 10° 

-0.056 

-0.175 

0.17 

0.011 

+0.006 


4 x 10° 

+0 . 726 

+0.636 

0.13 

0.0072 

0.773 


5 x 10° 

1.51 

1.44 

0.11 

0.0066 

1.56 


6 x 10° 

2.30 

2.24 

0.092 

0.0056 

2.34 


7 x 10° 

3.09 

3.04 

0.080 

0.0048 

3.18 


1 x ,10* 

5.48 

5.44 

0.057 

0.0033 

5.50 

- 

3 x 10 1 

21.40 

21.39 

0.020 

0.0011 

21.41 


7 x 10 1 

53.28 

53.28 

0.0085 

0.00048 

53.29 


1 x 10 2 

77.20 

77.20 

0.0060 

0,00033 

77.21 




Table VI. - Annulus-Groove Model 


Net -Leakage Solution x 10 


-14 


for Seal No. 


2 at Zero rpm 



N/AP 


0 

N/AP Composite Solutions 1 


1/N., No 

is. 

o Part. 
Corr , 

Part. Corr. 
1. 

Annulus 

Diffusion 

2 

Groove 

Diffusion 

3 

Composite Flow 
1 + 2 + 3 

-9 

1 x 10 

1.81 

0.06 

1.92 

0.72 

2.70 

3 x 10" 2 

1.82 

0 , 16 

1.69 

0,49 

2.34 

-9 . 

7 x -10 

1.84 

0.34 

1.41 

0.31 

2.06 

1 x 10" 1 

1.86 

0.48 

1.30 

0.24 

2.02 - 

-l 

3 x 10 

1.96 

0.95 

0.87 

0.11 

1,93 

7 x 10 _1 

2.16 - 

1.50 

0.57 

0.05 

2.12 

1 x 10° 

2.30 

1.79 

0.45 

0.035 

2.28 

3 x 10° 

3.28 

3.06 

0.19 

0.012 

3.26 

7 x 10 u 

5.22 

5.12 

0.09 

0.005 

5.22 

1 x 10 1 

6.68 

6 . 60 

0.064 

0.0036 

6.67 

3 x 10 1 

16.39 

16.36 

0.022 

0.0012 

16.38 

7 x 10 1 ’ 

35.81 

35.80 

0.0094 

0.0005 

35.81 

1 x 10 2 

50.37 

50.36 

0.0065 

0,00036 

50.37 



Table VII,. Annulus Groove Model 


Net Leakage 

Solutions 

x 10 ^ for Seal No. -2 

at 2500 rpm 




N/AP 


N/AP Composite Solution 

- 

k 

1/N No Part. 

Part., Corr. Annulus Diffusion 

Groove Diffusion 

.Composite Flow 


Corr. 

1 

2 

3 

1 + 2 + 3 

-2 






1 x 10 

-0.26 

-3.06 

1.75 

0.71 

-0.60 

3 x 10~ 2 

-0.32 

-2.84 

1.52 

0.47 

-0.85 

7 x 10“ 2 

-0.41 

-2.52 

1.28 

0.29 

-0.95 

1 x 10" 1 

-0.47 

-2.36 

1.16 

0.23 

-0.97 

3 x 10 _1 

-0.67 

-1.82 

0.79 

0.10 

-0.93 

7 x 10 -1 

-0.74 

5-1. 41 • 

0.51 

0.047 

-0.85 

1 x 10° 

-0.70 

-1 .<21 

0.41 

0.033 

-0.77 

2 x 10° 

-0.38 

-0.66 

0.23 

0.016 

-0.41 

3 x 10° 

+0.04 

-0.16 

0.17 

0,011 

+0.02 

4 x '10° 

0.48 

+0.32 

0.13 

0.0072 

0.46 

5 x 10° 

0.94 

0.81 

0.11 

0.0066 

0.93 

6 x 10^ 

1.40 

1.30 

0.092 

0.0056 

1.40 

7 x 10° 

1.86 

1.77 

0.080 

0.0048 

1.85 

1 x - 10 1 - 

3.27 

3.21 

0.057 

0.0033 

3.27 

3 x 10 1 

12.76 

12.73 

0.020 

0.0011 

12.75 

7 x 10 1 

31.78 

31.77 

0.0085 

0,00048 

31.78 

1 x 10 2 

46.04 

46.05 

0.0060 

0.00033 

46.06 
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4 ,2 Sealing Coefficient 

The sealing coefficient as defined in Equation (95) is a common 

dimensionless viscoseal performance index. The relationship which 

exists betwe'en 3ihe net leakage flow and the sealing coefficient can 

best be .shown by reference to Figure 16. A family of -net leakage 

curves as a function of the seal pressure ratio for a fixed shaft speed 

are shown,. The curves progress to the left for increasing values of 

pressure ratio until the (r + l)/(r - 1) factor approaches unity. 

p p 

When the family of net leakage flow curves are extended to the zero 
net flow condition, a corresponding sealing coefficient condition may 
be indicated for each of the pressure ratio curves . Operating the 
rarefied viscoseal at '’negative" net 'leakage is simply another way of 
saying that the viscoseal has now become a positive flow pump. 

. In Figures 17 through 22 sealing coefficient data and theoretical 
values are presented. Figures 17, 18, and 19 show results for 
Viscoseal No. 1 at 5000; 10,000; and 30,000 rpm. Figure 17 presents a 
comparison of experimental data and theoretical values as predicted 
by the Modified Hodgson Model. Figure 18 shoitfs the same type of 
comparison for the Slip Modified Boon and.Tal Model. Theoretical, 
values are shown for both the .mean-value and variable mean free path 
solutions .' Figure 19 is a similar comparison for the Annulus-, Groove 
Model., 

Figures 20 through 22 show experimental and theoretical values for 
Viscoseal No. 2. Again data and analytical values are shown for each 
of the three models -at .speeds of 5000; 10,000; and 30,000 rpm. 
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5. CONCLUSIONS 

Examination of Figures 10, 11, and 12.,' for Seal No. 1, indicates 
that the Annulus -Groove Model gives the best agreement with the 
experimental data - . The Slip Modified Boon and Tal Model gives the 
poorest agreement. All three models predict the general characteristic; 
for this low aspect (groove width ,to groove depth ratio) ratio, 1.002, 
seal. It is concluded that none of the theoretical models are completely 
adequate .for predicting the leakage through low aspect ratio viscoseals. 
The Annulus -Groove Model is recommended for low aspect ratio viscoseals 
and Inverse Knudsen Numbers greater than 5. 

For Seal No. 2, the results shown in Figures -13, 14, and 15 
indicate 'that the Slip Modified Boon and Tal Model gives a better 
correlation with experimental data than either of the other two models. 
The only exception is for static (zero rpm), conditions having Knudsen , 
Numbers greater than 0.5 where the Modified Hodgson Model shows better 
agreement. The aspect ratio of this geometry was 11.071 as .compared to 
1.002 for Seal No. -1, As -anticipated the Slip Modified Boon and Tal 
Model is much better for the higher aspect geometry than for low aspect 
ratio geometries. The variable mean free path solution is -only slightly 
better than the mean-value solution to the Slip Modified Boon and Tal 
Model. In view of .the added -complexity of the variable mean free path 
solution, it ;is .recommended that the mean-value Slip Modified Boon 
and Tal Model be used to predict viscoseal leakage for high 
aspect ratio geometries provided the^ Inverse Knudsen Numbers are 
greater than 5. It is of .interest .to note 'that the 2500 rpm data for 
Viscoseal No. 2 seem to approach a constant specific leakage value with 
increasing Knudsen number. 
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In Tables II through, VII a comparison between the continuum-slip 
and composite type -solutions are shown for the Slip Modified Boon and 
Tal Model and the Annulus-Groove Model. These particular cases were 
chosen for presentation .since they are the ones which gave maximum 
deviation between the two solution techniques. .Examination of Table'll 
shox'fs that the diffusion flow contributions in both the annular space 
and the sixteen, grooves are only one percent of 'the total composite 
flow at an Inverse Knudsen Number of 7.0 and decrease rapidly with 
decreased rarefi cation. At an Inverse. Knudsen Number of 0.01 it, can.be 
observed that the diffusion contributes 92 percent of .the total composite 
solution and thus .the transport is essentially free molecular. This 
same trend is 'evident for fhe other two static cases as shown in 
Tables IV and VI, In Tables III„ V, and VII- the .negative values .of 
specific leakage simply indicate that the viscoseal is acting as ,a 
pump, i.e., the viscous pumping flow “due to the grooves is .greater than 
leakage flow, ‘ 

For the 'Slip Modified Boon and Tal Model the maximum deviation 

between the composite solution and the .continuum solution is less than 

three percent for Inverse Knudsen Numbers gre.ater than 5.0. The deviation 

increases with' rarefication and has a maximum value of approximately 

twenty-five percent. This maximum occurs at 1/N - 0.01 and 1 zero rpm 

K 

for Seal No. 2. 

For the Annulus -Groove Model the maximum deviation is less .than. 

two percent for Inverse Knudsen Numbers greater than five. The maximum 

deviation was . approximately sixty percent" at 1/N “.0.01 and 2500 rpm 

K 


for Seal No, 2. 
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The good agreement between the rather complex composite solutions 
and the continuum-slip solutions for Inverse Knudsen Number greater than 
five is apparent. It is concluded that the complex composite solutions 
cannot be justified for Inverse Knudsen Numbers greater than 5.0. 

Consideration of the results shown in Figures 17 through 22 
reveals several .important aspects of the three analytical .models . Only 
the Modified Hodgson Model and the .variable mean free path solution -to 
the Slip Modified Boon and Tal Model predict the general sealing 
characteristics which are evidenced by the experimental data. In 
general the models predict better sealing performance (smaller sealing 
coefficients) than what actually existed. It may.be observed that the 
experimental data indicate that the ’Sealing coefficient parameter 
becomes speed sensitive as the sealant gas becomes rarefied. This 
sensitivity to, shaft, speed for the rarefied-gas viscoseal is unique 
from the -operation in the laminar continuum regime where the sealing 
coefficient is independent of shaft 'speed (14, 15). This continuum, 
characteristic may be seen by noting the asymptotic behavior of the 
experimental data for 5000 and 10,000 rpm as the continuum (1/N > 

100) regime is approached. These data may be compared since the mean 
radial clearance is essentially the same for both speeds. 

Examination of .Figures 18 and 21 reveals that the mean-value 
solution to the Slip Modified Boon and Tal Model can predict a seal 
pressure -differential, AP, which may exceed the magnitude of the upstream 
pressure. This solution would imply a .negative absolute value for the 
downstream pressure, a condition which cannot occur in reality. • A 
limiting criteria to be applied to this theoretical solution to avoid 
this condition is predictable by considering the variables involved. 
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As stated previously, the Knudsen number 


N K 


A 

c 


li(2TT^ T) l/2 
2Pc 


is based on the average pressure in the seal 


p + P 
r a T b 


2P t - AP 


As the downstream pressure approaches zero, then. AP approaches the 
upstream pressure and for this limiting condition 


— Ap 

p Hj 

r 't 2 


Thus, the Inverse Knudsen Number variable may be written ,as 


(i/n k ) - 


APc 


U(27T^'T) 1/2 


The Sealing Coefficient variable is defined as 


A - 6UUL 

** A O 

APc 

Solving for Ap in the Inyerse Knudsen Number relationship and 
substituting this into the above equation gives the limiting functional 
relationship 


A limit *' 


6UL 


R 


o „.l/2 
c(2tt— T) 


(1/n k ) 


(105) 


For a given seal geometry, gas, temperature, and operating speed, the 
limiting seal coefficient function varies .inversely with the. Inverse 
Knudsen Number, l/N^,. On a logarithmic plot such as Figure 17, this 
relation will ,plot as a straight line with a slope of minus unity. 
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Points which fall above and to the right of this line are possible while 
those which fall below and to the left are not possible in that, a 
negative absolute pressure downstream is implied. The limiting 
relation of Equation (105) is speed sensitive and thus the limit 
condition can be established for each selected shaft speed. Limit curves 
for speeds of 5000, 10,000, and 30,000 rpm from Equation (105) are 
shown on Figures 18, 19, 21, > and 22. The experimental .sealing 
coefficient data agree well with these limit curves and all points 
fall within the "possible" regions. In addition, the speed sensitivity 
of the experimental data agrees with that indicated by Equation (105) . 

It should be noted that the limiting function of Equation (105) 
results from a limit condition of the variables involved and applies 
to any and all viscoseal analyses using these variables. 

The limiting relationship just established is shown on Figures 18 
and 21 to indicate that the mean-value solution to the Slip Modified 
Boon and Tal Model x<?hen used in con j traction with Equation (105) is a 
useful technique for predicting sealing coefficient .performance. This 
limiting relationship can also be used with the Annulus-Groove Model. 

The devaiation of predicted performance from experimental results 
is approximately the .same for each of the three analytical models 
when the limiting relationship is used on the mean-value Slip Modified 
Boon and Tal Model and the Annulus-Groove Model. Theoretical values 
are not presented for the composite type solutions since there. is 
negligible difference between the composite solutions for sealing 
coefficient and the continuum-slip solutions. 

Since most practical seal geometries will result in a high 
aspect ratio, near 10, it is concluded that the simplicity of the mean- 
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value solution to the Slip Modified Boon and Tal Model in sconj unction 
with the limiting relationship justifies its use. For low aspect ratio 
geometries the Annulus-Groove Model in conjunction with the limiting 
relationship is recommended. The complexity of the variable mean free 
path solution (Slip Modified Boon and Tal Model) cannot be justified 
in view of the small improvement in- predicted , value.- The composite 
solutions are not practical for the .same reasons. 

6. RECOMMENDATIONS 

In general it is recommended that- the mean-value solution to the' 
Slip Modified Boon' and Tal Model, Equation (55), be used to predict 
viscoseal leakage for rarefied conditions resulting in Inverse Knudsen 
Numbers greater than 5. No satisfactory analytical model has been 
developed in this study to predict leakage for Inverse Knudsen Numbers 
less ,than> 5. 

To predict sealing performance it is recommended, in general, 
that the .mean-value solution to the Slip Modified Boon and Tal Model, 
Equation (50) , in -conjunction with the limiting relationship. 

Equation (105) , be used. 

It should be noted that the variable .mean-free-path solution to 
the Slip Modified Boon and Tal Model is .more accurate for all conditions. 
If use is made of a high speed digital computer its complexity becomes 
less important and thus its use may be justified. 
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Figure 1. Basic Elements of the Viscoseal. 
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Figure 4. Developed Viscoseal Geometry. 
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Figure 7. Viscoseal Test Section. 
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Figure 8. Schematic Diagram o£ the -Experimental Test Apparatus. 
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Figure 12. Net Leakage Flow Versus Inverse Knudsen Number 
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Figure 13 r Net Leakage Flow Versus Inverse Knudsen Number 
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Figure 14. .Net Leakage -Flow Vers-us Inverse Knudsen Number 
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Figure 16. Relation Between the Net Leakage Flow and the Sealing Coefficients 
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Figure 17. Sealing Coefficient Versus Inverse Knudsen Number 
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Figure 18. Sealing Coefficient Versus Inverse Knudsen Number 
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Figure 19. Sealing Coefficient Versus Inverse JCnudsen Numb 
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APPENDIX A. CONSTANTS FOR HODGSON'S MODEL 


The following expressions for the coefficients, in Hodgson’s 


Model are taken from reference (9). 


_ NFdb 2 h 

120 (b + h) c 

bh 3 (b + w)K 3 


12lJTTd 


SK^bh^Cb + w) 
C = 3ird(b + h) 

n , _ Tfdc 3 Cb + w) 
12yw 


,kT .1/2 
l 2irm j 


E „ 8K 2° 2 W0> + ») kT 1/2 

3w 7W 

K 2 (b + h) 2 h 

C i- — 

kT 1 /2 

8K i bu( W [2K l b " K 3 (b + ' h )3 

2K l b 

C 2 = K 3 (b + h) C 1 


C 3 = 


8l V 2K 2 - 


°4 ' 2K 2 C 3 


K - In (a + (1 + a 2 ) 1/2 ) + a 2 la A . + O * ^' 1 . 

8a 2 a 


+ ~ (1 + a 3 - (1 + a 2 ) 3 / 2 


)] 


K„ = 1 - 0.63 ~ tanh 
3 b 2h 



93- 


Values of K 2 

w/c 0.1 0.2 0.4 0.8 1 2 3 4 5 10 

K 2 0.036 0.068 0,13 0.22 0'.26 0.40 0.52 0,60 0.67 0.94 

3 

When w/c > 10, K£ = ^ In — . 


The groove flow coefficients, B, C, C-^, and C 2 are obtained 
from a consideration of flow in a long rectangular tube. The continuum 
groove, flow coefficient, B, is obtained, from the continuum, solution .for 
Poiseuille flow in a long , rectangular tube. From a consideration of 
purely molecular flow, the free-molecule flow coefficient, C, can be 
derived. The basic - free-molecule flow equation used by Hodgson is 
given in Reference (32) and is written as 


Q « 


16 

3 


K r kT \ in 


0 2 dP 
S dZ ' 


(106) 


Equation (106) is attributed to Knudsen (32, p. 35), but the constant 
K , which is equal to K~, is obtained by Clausing (32, p. 40). 0 and 

S are the groove ctoss sectional area and perimeter, respectively. 

The two remaining groove flow constants, and C 2> are obtained in 
the same way Knudsen determined V and £ in Equation (1) . As was 
pointed out in Section 2.1, the determination of the ratio V/C requires 
knowledge of the continuum with slip solution. Hodgson does not take' 
an exact approach at this .point since he determines C^/C 2 from the 
continuum with slip solution for flow between infinite parallel flat 
plates rather than flow in a long rectangular tube, 

Hodgson obtains the land leakage flow coefficients, D, E, C^, 
and from a consideration of flow in a thin slit-like tube. In this 
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case the continuum coefficient, D, is obtained from the continuum 
solution for Poiseuille flow between parallel .flat .plates. The same 
basic model, Equation (106), for free-molecule flow is used to determine 
the free-molecule flow .coefficient, E. now becomes K^, which is 
also attributed to Clausing. 0 and S now represent the cross sectional 
area of the land leakage passageway- and its perimeter, respectively. 

The constants and are also determined using the same approach 
Knudsen used to determine V and £ in Equation (1) „ For the continuum 
with slip solution, Hodgson once again uses the continuum with slip 
solution for flow between infinite parallel flat .plates which in this 
case is the proper solution to use. 
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APPENDIX B. TABULATED EXPERIMENTAL REDUCED DATA 



Net Leakage Reduced Data 


Data Point No. P^, 




1K13K 

1 

16,800 

7230 

2 

19,150 

7705 

3 

21,400 

8155 

4 

22,600 

8430 

5 

17,000 

7465 

6 

12,150 

565 

7 

9,370 

538 

8 

8,100 

522 

9 

6,850 

512 

10 

5,450 

500 

11 

4,180 

495 

12 

4,200 

490 

13 

2,950 

487 

14 

1,870 

482 

15 

13,800 

1250 

16 

1,740 

255 

17 

1,700 

260 

18 

1,720 

275 

19 

752 

227 

20 

428 

132 

21 

318 

129 

22 

251 

131 ' 

23 

216 

131 

24 

4,000 

131 

25 

2,320 

172 

26 

1,460 

166 

2,7 

1,025 

168 

28 

600 

172 

29 

475 

131 
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Seal No. 1, Zero rpm ~ Argon Gas 


i/n k 

N/AP 
x 10 

24.68 

18.16 

27.57 

20.30 

30.32 

21.60 

31.87 

22.08 

25.09 

17.45 

13.05 

10.60 

10.15 

9.06 

8.84 

8.24 

7.54 

7.28 

6.10 

6 . 34 

4.79 

5.50 

4.80 

5.46 

3.52 

4.65 

2.41 

3.81 

15.40 

11.31 

2.10 

3.47 

2.01 

3.60 

2.05 

3.59 

1.00 

3.23 

0.57 

3.02 

0.46 

2.83 

0.39 

2.40 

0.36 

2.48 

4.23 

5.26 

2.55 

4.15 

1.66 

3.46 

1.22 

3.30 

0.77 

2.77 

0.62 

3.20 
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Net Leakage Reduced Data for Seal No. 1, Dynamic Speeds - Argon Gas 








Data 

Point, 

No. 

P T 

Uhg 

P b 

yhg 

l/N 

K 

N/AP 
x 10“ 14 

c 

in. x 10 3 


10,000 rpm, c x 10 3 - 3.87 in,, r^ = 146.4 


1 

9,380 

335 

8.83 

0.558 

3.81 

2 

6,410 

60 

6.03 

0.117 

3.80 

3 

13,200 

86 

12.9 

1.50 

3.90 

4 

16,700 

93 

15.9 

2.74 

3.87 

5 

28,700 

150 

27.5 

7.80 

3.87 

6 

32,500 

166 

30.9 

9.80 

3.87 

7 

8,050 

64 

7.78 

0.237 

3.91 

10 

31,000 

233 

30.7 

10.15 

4.01 


5,000 

rpm c x 

10 3 = 3.85 in., 

r = 176 .6 

n 



11 

16,800 

63 

16 . 6 

7.29 

4.01 

12 

11,200 

42 

10.6 

4.14 

3.84 

13 

7,900 

33 

7.18 

2.61 

3.69 

14 

6,150 

30 

5.88 

1.88 

3.88 

15 

3,700 

27 

3.51 

0.953 

3.84- 

16 

1,950 

26 

2.03 

0.540 

3.76 

17 

1,250 

29 

1.32 

0.290 

3.94 


30,000 

rpm, c x 

10 3 = 3.29 in. , 

r * 807 
P 


18 

44,000 

53 

36.0 

0.451 

3.33 

19 

49,900 

64 

40.9 

1.'36 

3.33 

20 

39,800 

49 

31.4 

0.038 

3.20 
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Rarefied Sealing Coefficient Reduced Data for Seal No. 1 - Argon Gas 









Data Point 
No, 

P T ' 

P b 

_U h -g 

1/n k 

A 

Re 

c 

c 

in. x 10 3 



10 

,000 rpm. 

c x 10 3 

= 

3.98 ini 



SP1 

7,160 

267 

7.47 


40.2 

1.01 

4.18 

SP2 

9,860 

1220 

11; 1 


31. .8 

1.49 

4.15 

SP3 

5,200 

275 

5,01 


68.3 

0.67 

3.80 

SP9 

14,200 

3495 

16.9 


28.3 

2.28 

3.98 

SP10 

17,200 

6605 

22.6 


29.3 

3.02 

3.94 

SP13 

19,200 

7755 

25.5 


27.1 

3.42 

3.94 

SP14 

8,700 

44 

8.29 


36.4 

1.10 

3.91 

SP19 

7,350 

44 

6,96 


43.1 

0.93 

3.91 


5000 rpm, c 

x 10 3 = - 

3. 

99 in. 



SP4 

2,600 

166 

2.69 


60.8 

0.18 

4.04 

SP5 

11,000 

4730 

15.6 


22.5 

1.03 

4.11 

SP6 

3,880 

176 

3.94 


39.9 

0.26 

4.04 

SP7 

3,160 

407 

3.47 


53.7 

0.23 

4.04 

SP15 

2,850 

50 

2.72 


56.7 

0.18 

3.90 

SP16 

13,300 

8245 

20,6 


30.2 

1.39 

3.98 

SP17 

3,920 

146 

3.69 


46.5 

0.25 

3.73 

SP18 

4,400 

382 

4.65 


36 . 6 

0.31 

4.04 

SP20* 

4,200 

559 

4.47 


43.6 

0.30 

3.90 

SP21* 

7,000 

1835 

8.44 


29.6 

0,57 

3.97 

SP22* 

7,550 

1710 

8.94 


25.5 

0.60 

4.01 

SP23* 

13,100 

8145 

20.7 


29.8 

1.39 

4.04 

* Argon 

blanket 

on graphite ring 

seal. 




30,000 rpm, c x 10 3 

3 o 4l in o 



SP11 

35,200 

104 

29.3 


34.7 

11.8 

3.44 

SP12 

38,100 

68 

30.9 


33.5 

12,5 

3.37 



-99- 


Continuum Sealing Coefficient Tabulation Viscoseal No. 1, Gas-Air 


Data 

Point 

No. 

Speed 

rpm 

AP 

in.-H^O 

A 

R 

e 

c 

c 3 
in. x 10 



Test 

Date -6-25-69 



2 

4,905 

2.45 

23.7 

88.1 

4.18 

3 

8,080 

3.74 

25.6 

145.4 


4 

11,500 

5.15 

26.5 

207.3 


5 

13,200 

5.74 

27.2 

238.1 


6 

19,600 

8.49 

27.4 

354.7 


7 

23,400 

10.38 

26.6 

424.4 


8 

25,400 

11.45 

26.2 

461.3 


9 

29,600 

14.27 

24.5 

539.5 


10 

31,800 

16.50 

22.8 

581.2 


11 

34,600 

19.58 

20.9 

634.9 


12 

21,600 

9.45 

27.0 

391.4 


13 

9,700 

4.45 

25.8 

174.6 


15 

4,800 

2.40 

23.7 

87.1 


16 

6,400 

3.19 

24.5 

116.2 


17 

7,900 

3.69 

25.4 

143.5 


18 

9,700 

4.40 

26.0 

176.4 


19. 

II 56 OO 

5.12 

26.8 

211.1 


20 

13,150 

5.74 

27.1 

239.5 




Test 

Date 7-2-69 



22 

1,272 

0.73 

20.6 

23 . 1 

4.18 

23 

3,219 

1.69 

22.6 

58.5 


24 

4,644 

2.34 

23.5 

84.5 


25 

6,399 

3.08 

25.0 

116.5 


26 

7.869 

3.68 

25.2 

143.4 


27 

9,582 

4.35 

26.0 

174.7 


28 

11,562 

5.13 

26.7 

211.0 


29 

13,269 

5.83 

27.0 

242.5 


30 

7,392 

3.51 

24.9 

134.6 


31 

1,278 

0.74 

20.4 

23.2 
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Continuum Sealing Coefficient Tabulation Viscoseal No. 1, Gas 7 Air 



Data 

Point 

No. 

Speed 

rptti 

AP 

in.-H 2 0 

A 

R 

e 

c 

C 3 

in. x 10 


Test Date --9-2-69 - 


32 

1,566 

0.87 

21.5 

28.1 

33 

3,540 

1.80 

23.3 

63.7 

34 

5,130 

2.57 

23.6 

92.4 

35 

6,690 

3.17 

25.0 

120.6 

36 

7,980 

3.70 

25.5 

144.0 

37 

9,750 

4. 40 

26.3 

176.0 

38 

9,990 

4.50 

26.3 

180.4 

39 

13,800 

6.01 

27.2 

249.7 

41 

14,100 

6.21 

26.9 

255.2 

42 

19,230 

8.66 

26.3 

349.1 

43 

21,180 

9.34 

26.8 

384.8 

44 

23,400 

10.41 

26.6 

425.7 

45 

5,040 

2.50 

23.9 

90.8 

46 

1,278 

0.71 

21.3 

23.0 

47 

3,540 

1.81 

23.2 

63.7 

48 

5,010 

2.46 

24.1 

90.2 

49 

6,570 

,3.12 

25.0 

118.4 

50 

7,740 

3.60 

25.5 

139.6 

51 

9,720 

4.35 

26.4 

175.5 

52 

11,490 

5.04 

27.0 

207.6 

53 

13,830 

5.99 

27.4 

250.2 

54 

2,580 

1.34 

22.8 

46.4 

55 

510 

0 . 28 

21.6 

9.2 
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Continuum Sealing Coefficient Tabulation Viscoseal No. 1, Gas-Air 


Data 


ap 

A 

R 

Point, 

Speed 

in. H 2 0 


0 

c 

No. 

rpm 

.. 





Test 

Date 9-17-69 


66 

507 

0.28 

21.5 

9.2 

67 

2,460 

1.29 

22.6 

44.6 

68 

3,600 

1.84 

23.2 

65.3 

69 

5,010 

2.46 

24.1 

90.9 

70 

6,360 

3.04 

24.8 

115.5 

71 

7,470 

3.39 

26.1 

135.8 

72 

699 

0.39 

21.2 

12.7- 

73 

1,194 

0.66 

21.4 

21.6 

74 

1,338 

0.74 

21.4 

24.2 

75 

1,530 

0.85 

21.3 

27.7 

76 

1,710 

0.94 

21.5 

31.0 

77 

1,950 

1.05 

22. 0 

35.3 

78 

2,190 

1.17 • 

22.1 

40.0 

79 

270 

0.15 

21.3 

4.9 

80 

627 

0.35 

21.2 

11.4 

81 

867 

0.48 

21; 4 

15.7 

82 

762 

0,42 

21.5 

13.8 

83 

1,374 

0,75 

21.7 

24.9 

84 

1,647 

0.90 

21.7 

29.8 

85 

2,229 

1.19 

22.2 

40.4 

86 

2,970 

1.54 

22.8 

53.9 

87 

3,960 

2.00 

23.5 

71.8 

88 

4,980 

2.45 

24.1 

90.4 

89 

5,760 

2.78 

24.6 

104.6 

90 

7,140 

3.35 

25.2 

129.8 

91 

8,280 

3.80 

25.8 

150.6 

92 

9,090 

4.12 

26.2 

165.4 

93 

10,410 

4,63 

26.6 

189.5 

94 

11,640 

5.12 

26.9 

212.0 

95 

12,720 

5.57 

27.1 

231.8 

96 

13,530 

5.95 

27.0 

246.7 



- 102 - 


Net Leakage Reduced Data for Seal No. 2, Zero rpm, Argon Gas 



Data Point -No. 

yhg yhj 

t 

3 

N/AP 
x 10“ 14 


1 

17,000 

425 

18.00 

17.50 

2 

14,300 

403 

15.10 

15.40 

3 

10,200 

385 

10.90 

11. 80 

4 

8,080 

382 

8.68 

9.87 

5 

3,550 

375 

4.02 

6.12 

6 

1,500 

363 

1.91 

4.43 

7 

770 

360 

1.16 

4.05 

8 

2,460 

355 

2.89 

5.26 

9 

5,800 

355 

6.31 

8.13 

10 

440 

3.0 

0.45 

3.81 

11 

145 

1.8 

0.15 

3.68 

12 

251 

2.2 

0.26 

3.81 

13 

680 

4.0 

0.70 

3.92 

14 

65,500 

920 

68.10 

6.76 

15 

40,500 

420 

42.00 

3.97 

16 

28,900 

265 

29.90 

2.93 

17 

64 

1.6 

0.067 

3.58 
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Net Leakage Reduced Data for Seal No. -2, Dynamic Speeds, Arg'on Gas 


Data 






Point 

No. 

P T 

_ JM 

P b 

Dkg 

i/» K 

N/AP 

x-10 

c - 

3 

in. x 10 


2500 rpm, c x 10- ~^- v '31 9'4''‘ir£7', "r p ’ =' 21275' 


1 

21,200 

140 

20.70 

13.90 

3.95 

2 

15,700 

101 

15.20 

10.00 

3.' 91 

3 

11,000 

83 

10.40 

6,88 

3.83 

4 

9,000 

77 

8.62 

4.81 

3.87 

5 

6,200 

70 

5.95 

3.10 

3.87 

6 

4,670 

68 

4.68 

2.35 

4.03 

7 

3,240 

66 

3.27 

1.55 

4.03 

8 

1,850 

65 

1.89 

0.888 

4.03 

9 

29,200 

235 

28.50 

21.10 

3.95 

10 

42,900 

395 

42.40 

33.’ 90 

3.99 

11 

73,000 

978 

72.40 

63.70 

3.99 

12 

1,180 

2.1 

1.11 

1.14 

3.83 

13 

222 

0.6 

0.20 

1..02 

3.71 

14 

97 

0.45 

0.093 

1.017 

3.90 

15 

59 

0.40 

0.055 

1.014 

3.79 

16 

12,450 

53 

12.10 

7.63 

3.95 

17 

2,840 

6.0 

2 .67 

1.77 

3.83 

18 

1,410 

2.4 

1.33 

1.15 

3.83 

19 

952 

1.7 

0.92 

1.09 

3.95 

20 

505 

1.0 

0.48 

1.01 

3.83 


30,000 

_ 3 

rpm, c x 10 = 

3.51 in. , r 
P 

= 246.8 


21 

132,500 

1050 

120.5 

35.7 

3.68 

22 

103,700 

530 

94.1 

16 . 1 

3.68 

23 

74,000 

150 

62.6 

1.73 

3.44 

24 

88,500 

515 

70.3 

5.19 

3.22 
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Net Leakage Reduced Data for Seal No. 2, Dynamic Speeds, Argon Gas 



Data 

Point 

No. 

P T 

ykg 

P b 

Dbg 

!/N k 

N/AP 
x 10 

c 

in. x 10 3 


5000 rpm, c x 10 3 - 3.92 in., *» 891.9 


25 

29,800 

200 

29.40 

15.6 

3.-99 

26 

44,500 

360 

43.50 

28.2 

3. '95 

27 

78,800 

1000 

77.30 

59.0 

3.95 

28 

4,100 

2.9 

3.61 

0.509 

3.59 

29 

2,850 

1.5 

2.59 

0.352 

3.71 

30 

6,660 

8.5 

6.20 

1.25 

3.79 

31 

1,670 

0.75 

1.45 

0.199 

3.54 

32 

525 

0.25 

0.41 

0.111 

3.19, 

33 

760 

0.30 

0.65 

0.199 

3.47 

34 

388 

0.20 

0.34 

0.157 

3.59 

35 

355 

0.20 

0.31 

0.163 

3.54 

36 

25,100 

137 

24.70 

11.7 

3.99 

37 

15 , 800 

81 

15.70 

5.96 

4.03 

38 

9,500 

59 

9.26 

2.00 

3.95 

39 

19,100 

106 

18 . 10 

8.1'4 

3,97 

40 

12,000 

80 

11.40 

3.64 

3.77 

41 

6,330 

69 

6.30 

0.791 

3.77 

42 

5,200 

75 

4.97 

0.561 

3.75 


10 

,000 rpm, c x 

10 3 = 3.85 

in. , r = 1879.2 
P 


43 

17,100 

70 

16.3 

1.11 

3.87 

44 

27,800 

109 

26.8 

5.72 

3.91 

45 

33,500 

135 

31.9 

'8.55 

3.86 

46 

23,100 

115 

21.8 

3.23 

3.83 

47 

7,150 

0.6 

6.61 

0.0251 

3.77 

48 

14,600 

5.5 

12.40 

0.331 

3.46 

49 

16 ,500’ 

11.0 

14.20 

0.648 

3.51 

50 

10,900 

2.2 

9.71 

0.147 

3.63 

51 

88,100 

1000 

87.10 

52.7 

3.99 

52 

44,400 

265 

42.80 

16.5 

3.91 

53 

36,100 

180 

34.80 

10.7 

3.91 

54 

66,000 

550 

64.50 

33.0 

3.95 
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Rarefied Sealing Coefficient Reduced Data for Seal No. 2' t Argon Gas 


Data 

Point 

No. 

p 

T 

Uhg 

P b 

Uhg 

1/h k 

A 

R 

e 

C f 

c 

in. x 10 3 


5000 

rpm, c-x 10. 

■= 3.92 

1 



SP1 

4,600 

90 

4.41 

35.0 

0.296 

3.91 

SP2 

35,000 

28,000 

59.80 

22.2 

4.006 

3.94 • 

SP3 

33,000 

25,800 

55.77 

21.6 

3.739 

3.94 

SP4 

12,000 

4,400 

15.70 

20.0 

1.053 

3.98 . 

SP5 

24,000. 

16,600 

38.90 

20.6 

2.608 

3.98 

SP6 

8,950 

2,330 

10.80 

23.0 

7.245 

3.98, 

SP7 

8,250 

1,670 

9.24 

24.4 

6.195' 

3.87 

SP8 

3,180 

50 

2.95 

53.7 

0.-198 

3.79 

SP14 

42,000 

34,800 

73.80 

21.0 

4.945 

3.99 

SP15 

56,200 

49,000 

102.80 

20.3 

6.893 

4.06 

SP18 

2,610 

66 

2.36 

70.7 

0.158 

3.66 

SP19 

520 

3.8 

0.45 

364 

0.030 

3.58 


10,000 rpm, c x 10 3 = 3.1 

36 in. 



SP9 

39,000 

23,000 

58.20 

19.7 

7.804 

3.90 

SP10 

14,000 

695 

13.54 

24.6 

1.’817 

3.83 

SP11 • 

10,700 

51 

9.91 

30.8 

1.329 

3.83 

SP12 

19,000 

3,380 

20.80 

20.7 

2.788 

3.86 

SP13 

24,900 

8,160 

31.00 

18.9 

4.162 

3.90 

SP16 

54,000 

36,200 

87.30 

16.8 

11.704 

4.02 

SP22 

4,450 

3 

3.87 

83.1 

0.519 

3.61 


30, 

000 rpm, c x 

10 3 = 3. 

.23 in. 



SP20 

111,200 

42,200 

123.00 

18.8 

49.463 

3.33 ■ 

SP21 

50,200 

1.1 

36.10 

32.1 

14.534 

2.99 

SP17 

63,000- 

205 

51.1 

20.3 

20.564 

3.36 
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Continuum Sealing Coefficient Tabulation Viscoseal No. 2, Gas-Air 



Data 

Point 

No. 

Speed 

rpm 

AP 

in. H 2 0 

A 

R 

e 

c 

c 

in. x 10^ 


1 

543 

0.48 

13.43 

9.73 

4.18 

2 

1,148 

0.99 

13.77 

20.6 

4.18 

3 

4,140 

3.50 

14.31 

73.7 

4.14 

4 

3,480 

2.98 

14.,13 

61.9 

4.14 

5 

2,850 

2.48 

13.91 

50.7 

4.14 

6 

2,310 

2.00 

13.98 

41.4 

4.14 

7 

1,710 

1.50 

13.80 

30.4 

4.14 

8 

1,110 

0.99 

13.57 

19.7 

4.14 

9 

570 

0.52 

13.27 

10.11 

4.14 

10 

4,800 

4.01 

14.49 

85.5 

4.14 

11 

5,430 

4.49 

14.70 

96.6 

4.13 

12 

6,090 

5.00 

14.88 

108.1 

4.12 

13 

6,690 

5.46 

15.04 

118.6 

4.11 

14 

7,380 

6.01 

15; 15 

130.6 

4.10 

15 

8,040 

6.51 

15.31 

142.0 

4.09 

16 

8,700 

7.02 

15.51 

153.0 

4.07 

17 

9,300 

7.47 

15.90 

162 .'0 

4.03 

18 

10,020 

8.01 

16.21 

173.4 

4.00 

19 

10,590 

8.48 

16.51 ■ 

181.5 

3.96 

20 

11,250 

9.01 

16 .'52 

193.0 

3.96 

21 

11,790 

9.46 

16.48 

202.4 

3.96 

22 

12,510 

10.'05 

16.46 

214.9 

3.96 

23 

13,110 

10.54 

16.45 

225.3 

3.96 

24 

14,520 

11.74- 

15.80 

254,3 

4.03 

25 

13,680 

11.07 

15.78 

239.4 

4.03 





-107- 


Continuum Sealing Coefficient Tabulation Viscoseal No, 2, GasrAir . 


Data 

Point Speed AP A R c 

No. ' rpm in.-H^O e c in. x Hr 


Test Date-- 9-2-70 


26 

4,950 

4. 11 

14.30 

89.1 

4.18 

27 

10,200 

8.00 

15,. 28 

183.6 

4.16 

28 

14,850 

11.8-3 

15.19 

267.2 

4.14 

29> 

16,440 

13.31 

14.95 

296.4 

4.14 

30 

18,180 

15.01 

14.66 

328.5 

4.14 

31 

19,590 

16.54 

14.33 

354.6 

4.14 

32 

21,210 

18.34 

14.00 

384.8 

4.14 

33 

22,650 

20.14 

13.61 

411.9 

4.14 

34 

23,700 

21.07 

13.24 - 

431.7 

4.14 

35 

24,930 

23.26 

12.97 

455.0 

4.14 

36 

25 ,920 

24.79 

12.65 

474.0 

4.14 

37 

27,000 

26.55 

12.31 - 

494.8 

4.14 

38 

27,900 

28.08 

12.02 

512.3 

4.14 

39 

213,890 

29.94 

11.66 

531.-7 

4.14 

40 

29,760 

31.81 

11.32 

549.0 

4.14 

41 

30,510 

33.41 

11.05 

563.9 

4.14 

42 

31,410 

35.58 

10.67 

582.1 

4.14 

43 

32,130 

37.32 

10.41 

596.7 

4.14 

44 

32,940 

39.83 

10 .12. 

613.2 

4.14 

45 

34,110 

42.43 

9.72 

637.3 

4.14 




